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ABSTRACT 

The  paper  studies  Loo(IR'^)-norm  approximations  from  a  space  spanned  by  a  discrete  set  of 
translates  of  a  basis  function  4>.  Attention  here  is  restricted  to  functions  (f>  whose  Fourier  transform 
is  smooth  on  1R'^\0,  and  has  a  singularity  at  the  origin.  Examples  of  such  basis  functions  are  the 
thin-plate  splines  and  the  multiquadrics,  as  well  as  other  types  of  radial  basis  functions  that  are 
employed  in  Approximation  Theory.  The  above  approximation  problem  is  well-understood  in  case 
the  set  of  points  E  used  for  translating  (j)  forms  a  lattice  in  IR*^,  and  many  optimal  and  quasi-optimal 
approximation  schemes  can  already  be  found  in  the  literature.  In  contrast,  only  few,  mostly  specific, 
results  arc  known  for  a  set  E  of  scattered  points. 

The  main  objective  of  this  paper  is  to  provide  a  general  tool  for  extending  approximation 
schemes  that  use  integer  translates  of  a  basis  function  to  the  non-uniform  case.  We  introduce  a 
single,  relatively  simple,  conversion  method  that  preserves  the  approximation  orders  provided  by 
a  large  number  of  schemes  presently  in  the  literature  (more  precisely,  to  almost  all  “stationary 
schemes").  In  anticipation  of  future  introduction  of  new  scheme  for  uniform  grids,  an  effort  is 
made  to  impose  only  a  few  mild  conditions  on  the  function  (f>,  which  still  allow  for  a  unified  error 
analysis  to  hold.  In  the  course  of  the  discussion  here,  the  recent  results  of  [BuDL]  on  scattered 
center  approximation  are  reproduced  and  improved  upon. 
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Radial  basis  function  approximation: 
from  gridded  centers  to  scattered  centers 

Nira  Dyn  k  Amos  Ron 


1.  Introduction 
1.1.  General 

The  terminology  “Radial  Basis  Function  Approximation”  usually  refers  to  approximation  to 
data  defined  on  a  set  E  of  scattered  points  in  IR**,  d  >  1  (referred  to  as  “centers”)  by  linear 
combinations  of  a  discrete  set  of  translates  (whose  standard  choice  is  the  original  E)  of  “a  radial 
basis  function"  <p.  Typical  examples  of  a  basis  function  (j)  are  the  bivariate  thin-plate  spline 

X  |x|^  log  lx|, 


and  the  Hardy's  muUiquadric 

X  I-+  (lx|^  + 

(where  here  and  hereafter  jxj  stands  for  the  Euclidean  norm  of  the  vector  x  €  IR**).  The  above 
functions  are,  indeed,  radially  symmetric  (i.e.,  obtained  by  composing  a  univariate  function  with  the 
map  X  |x|),  a  feature  which  greatly  facilitates  computations  with  such  functions.  However,  the 
title  “radial  basis  functions  approximation”  is  somewhat  misleading:  most  of  the  radially  symmetric 
functions  are  not  considered  as  candidates  for  the  basis  function  <f),  while,  on  the  other  hand,  certain 
classes  of  non-radially  symmetric  functions  do  fit  in.  We  elaborate  on  the  actual  conditions  required 
of  <t>  later  on  in  the  present  paper,  and  also  refer  the  reader  to  [R]  (particularly,  page  256  there). 

While,  as  alluded  to  before,  the  major  force  of  radial  function  approximation  is  in  the  area  of 
scattered  data  approximation  (and  particularly  interpolation,  cf.  the  surveys  [Dl,2],  [P],  [Bu4]), 
M.J.D.  Powell  and  I.R.H.  Jackson  (of  Cambridge,  England)  initiated  several  ye.*  s  ago  the  study  of 
approximation  from  a  principal  shift-invariant  space  (PSI,  for  short)  5  ^  erated  by  a  radial 
basis  function  <f>.  Such  a  space  5  is,  by  definition,  the  smallest  closed  (in  son.  Jerlying  topology) 
space  that  contains  all  shifts  (i.e.,  integer  translates)  of  the  basis  function  <p.  Together  with  the 
deepening  of  the  understanding  of  the  more  mathematically  accessible  uniform  case,  emerged  the 
hope  that  the  tools  and  observations  that  are  found  and  made  in  that  latter  case  will  serve  for  a 
better  understanding  of  the  practically  important  scattered  situation.  Indeed,  the  present  paper  is 
one  of  the  first  applications  of  the  now-established  theory  of  the  uniform  grid  case  to  the  scattered 
case. 

The  approximation  properties  of  PSI  and  other  shift-invariant  spaces  generated  by  radial  basis 
functions  are  analysed  in  detail  in  a  sequence  of  papers  written  by  Powell’s  students  I.R.H  Jackson 
and  M.D.  Buhmann,  by  the  two  of  us  and  by  several  other  authors  (cf.  the  above  mentioned  surveys 
and  the  references  within.  Specific  results,  hence  references,  which  are  relevant  to  this  paper  are 
sketched  in  §3).  Moreover,  those  results  on  approximation  orders  of  PSI  spaces  generated  by  radial 
basis  functions  served  as  one  of  the  major  forces  behind  the  general  results  on  approximation  orders 
of  [BR]  and  [BDeR].  In  fact,  the  totality  of  results,  techniques  and  observations  concerning  the 
approximation  properties  of  PSI  (and  related)  spaces  can  be  fairly  characterized  as  a  “solid  theory”, 
at  least  from  a  qualitative  point-of-view. 
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While  the  present  state-of-art  in  the  area  of  approximation  from  PSI  spaces  is  satisfactory,  the 
same  cannot  be  said  about  the  scattered  case,  (unless  the  spatial  dimension  is  1).  As  a  matter  of 
fact,  the  recent  paper  [BuDL]  contains  the  first  successful  application  of  tools  and  results  from  the 
shift-invariant  case  to  the  scattered  situation.  It  concerns  the  approximation  by  scattered  shifts  of 
any  one  of  the  basis  functions  of  [DJLR],  and  employs  quasi-interpolation  as  the  appro.ximation 
scheme.  Apparently,  the  theory  of  the  shift-invariant  case  is  invoked  in  that  paper  only  for  the  error 
analysis  part  (more  precisely,  in  the  polynomial  reproduction  argument),  and  not  for  the  (quite 
involved)  construction  of  the  approximation  map  there,  but,  our  post-analysis  of  the  map  from 
[BuDL]  reveals  that  that  construction  can  be  understood  as  a  specific  instance  of  the  following 
general  three-step  method.  We  use  here,  for  some  discrete  E  C  IR*^,  the  notation 


5h{^) 

to  denote  the  “span”  of  {<(>{•  —  4)}^6H1  where  the  precise  meaning  of  “span”  is  yet  to  be  defined. 
Step  1:  For  each  a  €  the  translate  —  a)  is  approximated  by  a  combination  (pa  of  {0(-  — 
^)}^gH-  This  gives  rise  to  a  linear  map  A  from  C~  into  (with  6^  the  sequence  which  is  1  at  ^ 
and  0  on  E\^,  A(^^)(a)  is  the  coefficient  A(^,a)  of  (f){-  —  in  the  definition  of  <pa-  A  then  extends 
by  linearity  to  the  extent  that  this  is  possible). 

Step  2;  An  approximation  scheme  L  that  is  used  to  approximate  from  S2idi<p)  is  modified  to  obtain 
an  approximation  map  La  into  S‘=i<p)  by  simply  replacing  each  appearance  of  (f>{-  —  a)  in  L  by  4>a- 
Step  3;  The  L  that  corresponds  to  the  specific  La  of  [BuDL]  (the  localization-quasi-interpolation 
map  of  [DJLR])  employs,  for  computing  the  approximant  Lf  of  /,  the  values  of  /  on  ZZ*^.  Since, 
presumably,  the  available  information  on  /  is  /|g,  the  map  A  (from  (a))  is  applied  to  /j,,  to  obtain 
approximations  {/aloga-'  to  {/(o()}a€a<‘-  The  actual  approximation  map  of  [BuDL]  could  then 
be  obtained  from  La  by  replacing  each  appearance  of  /(a)  by  fa- 

The  main  objective  of  this  paper  is  to  provide  a  general  tool  for  extending  approximation  maps 
that  use  uniform  translates  of  a  basis  function,  to  the  non-uniform  case.  We  do  that  by  showing 
that  Steps  1  and  2  above  can  be  combined  to  provide  a  general  method  for  converting  a  map  L  that 
acts  into  the  shift-invariant  S-^d{(}>)  to  a  map  La  that  approximates  from  (i.e.,  maps  into)  S-={4>) 
and  has  approximation  power  similar  to  that  of  L.  This,  quite  certainly,  cannot  be  achieved  without 
imposing  some  restrictions  on  the  underlying  basis  function  <f>:  ive  always  assume  throughout  the 
paper  that  the  function  (p,  considered  as  a  distribution  in  P'(]R‘*),  is  tempered,  and  that  its  Fourier 
transform  <p  coincides  on  IR“\0  with  some  continuous  function  (still  denoted  by  (p)  while  having  a 
certain  type  of  singularity  ( necessarily  of  finite  order)  at  the  origin.  In  turn,  such  an  assumption 
allows  us  to  construct  the  map  A  of  Step  1  independently  of  the  specific  <p  involved.  Indeed,  instead 
of  approximating  <p(-  —  a)  by  combinations  of  {^(-  —  0}fes>  will  actually  approximate  at  the 
origin  the  exponentials 

Cq  ;  X  e*“  *,  a  €  S'* 

by  combinations  of  under  the  requirement  that  A  is  bounded  as  a  map  from  ^oo(E)  to 

looiTZi'^)-  This  imposes  constraints  on  the  geometry  of  the  center  set  E  to  which  our  method  applies. 
The  problem  of  constructing  such  a  map  A  is  already  considered  in  [BuDL]. 

As  one  should  realize  from  the  above  description,  we  do  not  aim  at  establishing  approximation 
orders  for  spaces  generated  by  scattered  translates  of  a  suitable  basis  function,  nor  do  we  intend 
to  introduce  approximation  maps  especially  designed  for  dealing  with  scattered  approximation. 
Instead,  we  present  here  a  general  method  that  allows  us  to  convert  any  known  approximation 
scheme  on  uniform  grids  to  the  non-uniform  case,  while  preserving  (to  the  extent  that  this  is 
possible)  the  approximation  orders  known  in  the  former  case.  In  anticipation  of  future  introduction 
of  new  schemes  for  uniform  grids,  we  made  an  effort  to  impose  only  few  (mild)  assumptions  on 
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^  (in  addition  to  the  basic  one  specified  in  the  previous  paragraph).  Remarkably,  our  conversion 
method  based  on  Steps  1,2  above  applies  to  a  large  number  of  cases  presently  in  the  literature,  and 
our  error  analysis  is  also  quite  unified.  In  one  regard  we  sacrifice  generality  for  the  sake  of  brevity 
and  readability:  only  uniform-norm  approximations  are  considered  in  the  present  paper.  We  also 
mention  that  the  tools  developed  here  do  not  seem  to  apply  to  non-stationary  approximations  (such 
as  the  ones  that  employ  the  Gaussian  kernel,  etc.,  cf.  [MN],  [BuD],  [BR]  and  [BeLi])  that  lead  to 
spectral  approximation  orders. 


1.2.  The  setup  and  a  typical  result 

All  functions  considered  in  this  paper  are  either  real-  or  complex-valued  functions  defined  on 
IR**,  for  some  fixed  d.  Norms  of  such  functions  /  are  denote  by  ||/||...,  with  the  default  norm  in  this 
paper  is  that  of  Loo(lR‘^),  i-e., 

ll/ll  :=  ll/lloo  :=  II/IIl„(R-)- 

Norms  on  IR*^  itself  are  denoted  by 

\vl„,  V  €  IR'', 

with  the  default  norm  is  the  Euclidean  one: 

\v\^ 

Given  an  index  set  H,  and  a  function  set  we  set 

0s  :=  {0€h6S- 

In  the  special  case  when  E  C  IR*^,  and  the  function  set  is  the  set  of  all  translates  of  a  fixed  <f>  by 
^  €  E,  we  use  the  more  specific  notation  s(0): 

s(0):={0(--O:  4eH}. 

Given  a  function  set  F,  we  let 

SF 

be  the  closure  in  the  topology  of  uniform  convergence  on  compact  sets  of  the  algebraic  span  of  F. 
For  example,  this,  in  particular,  defines  precisely  the  space  53(0),  which  is  previously  introduced 
in  §1.1. 

We  are  interested  in  approximating  smooth  functions,  particularly  in  the  space 

k€2Z+ 

of  all  functions  whose  derivatives  of  orders  <  k  are  bounded.  Some  of  our  results  deal  with 
approximation  to  functions  in  the  larger  homogeneous  Sobolev  space 


of  all  functions  for  which  the  semi-norm 

|Q|,=fc 
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is  finite. 

As  is  almost  always  the  case  in  approximation  theory,  polynomials  and  exponentials  are  heavily 
employed  here.  We  have  already  mentioned  the  symbol 

C{,  4  €  IR**, 

for  the  complex  exponential  with  frequency  We  use  the  symbol 

n* 

to  denote  the  space  of  all  polynomials  of  degree  <  k. 

Finally,  given  /  continuous  on  1R‘^\0,  we  say  that  /  has  a  singularity  of  order  k'  at  the 
origin  if 

is  bounded  above  and  below  by  positive  constants  in  some  origin-neighborhood. 

For  a  given  H  C  IR'^,  our  interest  is  in  approximation  schemes  from  5h(<^)  that  are  obtained 
by  modifying  a  known  approximation  scheme  from  and  that  preserve  to  the  best  possible 

extent  the  approximation  power  of  that  latter  scheme. 

A  typical  approximation  scheme  considered  here  is  of  the  form 

La-/*-*  V'a A(/)(q'), 

where  the  functions  (^a)o€-Z‘‘  certain  functions  in  Sv{(f>),  known  to  satisfy  (at  least)  the  fol¬ 
lowing  “boundedness  condition” 

(1.2.1)  ^  e  Xoo(lR‘'). 

A  is  always  taken  to  be  a  bounded  operator  from  £oo(lR‘^)nC(IR‘^)  into  itself,  hence,  in  particular, 
the  sequence  ((A/)(a))£,g2''  bounded,  and,  hence  (due  to  (1.2.1))  La  is  well-defined  and  bounded 
from  Zoo(IR‘*)  n  (7(51*^)  into  S={(f>). 

The  actual  construction  of  La  starts  with  a  given  a  approximation  scheme  L  from  S2i,d{<f>)  of 
the  form 

L:f^  ^(.  -  q)  A(/)(a), 

with  A  as  above,  and  with  the  function  ^  being  a  linear  combination  of  the  shifts  of  the  original 
basis  function 

oez'' 

To  simplify  the  discussion  at  this  introductory  stage,  we  assume  that  the  sequence  ft  is  of  finite 
support  (the  actual  conditions  imposed  on  n  in  the  paper  are  much  milder).  At  a  minimum,  the 
function  ^  is  assumed  to  satisfy  the  condition 

(1.2.2)  Y.  IV'(  -«)|G^oo(]R‘'), 
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which  guarantees  the  map  L  to  be  well-defined  and  bounded  on  XooC®-**)  ^  C(SV^).  The  reader 
should  bear  in  mind  that  at  the  core  of  this  paper  is  the  assumption  that  the  map  L,  hence  the 
sequence  p,  the  function  •0,  and  the  operator  A,  are  all  given,  and  the  approximation  properties 
of  this  map  are  known,  too.  Our  goal  is  to  find  an  appropriate  selection  method  for  the  functions 
(^a)o<€Z‘‘  approximation  properties  similar  to  those  of  L.  Since  in  many  cases  the 

operator  A  is  the  identity,  we  simplify  further  our  introductory  discussion  here,  by  adopting  such 
an  assumption. 

Our  construction  of  the  functions  (ipa)a  is  done  as  follows.  We  first  approximate  each  <f>{--a), 
a  €  ZZ*^,  by  a  linear  combination 

(1.2.3)  K  “  53  Aia,  0<i>i-  -  ^)- 

€eH 

The  reader  who  feels  uneasy  about  the  lack  of  discussion  of  the  convergence  of  this  last  sum,  may 
assume  without  great  loss  that,  for  each  a  6  7Z‘^ ,  the  sequence  (A(a, is  finitely  supported. 
We  then  define  the  functions  (rpa)a  by 

=  53  4>0P{^-ot)- 

The  basic  criterion  of  the  approximation  properties  of  L  and  Lj\  that  we  employ  in  this  paper 
is  that  of  approximation  orders.  We  say  that  a  map  L  provides  an  approximation  order  A;  >  0 
for  a  smoothness  space  W,  if,  for  every  f  &  W, 

(1.2.4)  lkJ-I(cr„/)|l  =  0(/i*), 
with  Oh  the  scaling  operator 

f{h  ). 

Note  that,  since  the  uniform  norm  is  invariant  under  dilations,  we  could  write  (1.2.4)  equivalently 
as 

which  is  more  frequently  used  in  the  literature  for  defining  the  notion  of  approximation  orders. 

We  present  now  a  result  which  is  a  variant  of  Theorem  2.2.16,  and  which  can  be  regarded  as 
a  prototype  for  the  main  results  of  this  paper. 

Theorem  1.2.5.  Assume  that  the  basis  function  <p  is  tempered  (as  distribution)  and  that  its 
(distributional)  Fourier  transform  is  continuous  on  ]R‘*\0,  has  a  singularity  at  the  origin,  and  does 
not  vanish  identically  on  27rZ''\0.  Let  i>,  (<^a)a,  (V’a)a€Z'‘>  be  defined  as  above,  with 

A  =  identity.  Assume  that: 

(a)  p  is  Bnitely  supported. 

(b)  ^  satisfies  (1.2.2),  and  ^(0)  ^  0.  (Ip  is  continuous,  hence  well-defined  pointwise,  because  of 
(1.2.2).) 

(c) 

Then  (1.2.1)  holds  and  the  hence-well-defined  La  provides  (to  sufficiently  smooth  functions)  at 
least  the  same  approximation  order  as  provided  by  L. 

Here,  “sufficiently  smooth  functions”  are  functions  in  W^(1R‘^),  with  k  depending  on  the 
underlying  approximation  order. 
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The  outline  of  the  paper  is  as  follows.  §2  is  devoted  to  the  developments  of  the  general  theory. 
The  core  of  its  discussion  are  §2.1,  §2.2,  and  §2.7.  In  §2.1  we  impose  conditions  on  p,  ip  and  A  which 
are  sufficient  for  the  development  of  the  method  of  this  paper.  The  main  issue  in  that  subsection  is 
localization  which  amounts  to  the  proper  construction  of  the  function  ip  and  the  approximation  map 
L.  §2.2  discusses  the  extension  of  L  to  the  scattered  case,  i.e.,  the  construction  of  the  “pseudo-shifts” 
<pa  and  ipa,  and  the  introduction  of  our  approximation  map  It  also  includes  the  basic  error 
analysis  for  La  (cf.  Theorems  2.2.9  and  2.2.16).  In  §2.3,  we  modify  the  previous  error  analysis,  in  a 
way  that  relaxes  the  requirements  on  p  imposed  in  §2.2,  but  applies  to  a  smaller  set  of  approximands 
/.  Since  the  standard  results  in  the  literature  show  that  T  is  a  good  approximation  scheme  not 
only  for  smooth  bounded  functions,  but  also  to  functions  which  are  not  bounded,  but  some  of 
their  derivatives  are  bounded,  we  show  in  §2.4  that,  under  some  additional  assumptions,  such 
approximation  properties  are  inherited  by  La-  The  issue  of  polynomial  reproduction  is  at  the  heart 
of  error  analyses  which  are  based  on  the  quasi-interpolation  argument,  but  plays  no  real  role  in  the 
present  paper.  However,  we  show  in  §2.5  that  La  reproduces  polynomials  to  the  expected  degree, 
whenever  such  property  can  make  sense.  §2.8  suggests  some  relaxations  in  the  various  conditions 
assumed  in  §2.2  under  which  the  main  theorems  of  §2.2  are  still  valid.  Finally,  in  §2.7  the  following 
very  important  question  is  considered:  what  properties  of  the  matrix  A  imply  the  satisfaction  of 
the  crucial  condition  (c)  in  the  above  stated  theorem. 

Section  3  is  mainly  devoted  to  the  application  of  the  general  theory  to  specific  instances.  In 
§3.1,  some  of  the  approximation  schemes  and  approximation  order  results  of  [BR]  are  extended 
to  the  scattered  case.  The  quasi-interpolation  schemes  of  [DJLR]  are  extended  in  §3.3,  while  the 
extension  of  the  interpolation  schemes  of  [Bul,2)  is  the  issue  of  §3.2.  In  addition,  we  provide  in 
§3.4  a  detailed  discussion  of  the  scattered-center  approximation  schemes  of  [BuDL],  in  which  we 
approach  those  schemes  via  a  course  different  from  the  original  one  suggested  in  that  reference. 
In  this  way,  we  provide  an  alternative  error  analysis  of  the  [BuDL]-likc  schemes  that  is  based  on 
the  results  of  §2.2  here,  thereby  providing  extensions  for  the  [BuDL]  results  in  several  different 
directions. 

2.  The  approximation  maps  and  their  error  analysis 
2.1.  Localization  and  approximation  maps  on  uniform  grids 

A 

As  mentioned  before,  we  assume  that  the  Fourier  transform  0  of  the  basis  function  <p  is  regular 
and  even  continuous  on  1R‘*\0,  but  is  singular  at  the  origin,  and  such  assumption  implies  that 
<P{x)  cannot  decay  fast  to  zero  as  |x|  — »  oo;  for  instance  <p  0  7^1(11*^).  In  fact,  in  most  examples 
l</>(x)|  grows  polynomially  as  |il  tends  to  oo.  This  makes  it  hard  (although  not  impossible)  to  write 
explicit  approximation  schemes  which  are  based  on  shifts  of  <p.  The  standard  way  to  circumvent 
this  difficulty  is  via  a  localization  process,  which  uses  the  shifts  of  (p  to  construct  a  new  function  ip 
with  favorable  decay  properties  at  oo,  and  subsequently  implements  approximation  schemes  which 
are  b^lsed  on  the  shifts  of  ip. 

Formally,  the  localization  ip  can  be  written  as  a  linear  combination  of  shifts  of  <p  with  coefficients 
p  :  ZZ**  -»  C: 

(2.1.1)  ^“51 

aeZZ'* 

Usually,  the  localization  sequence  fi  decays  fast  at  oo  and  the  above  sum  converges  uniformly 
on  compact  sets  (as  a  matter  of  fact,  p  is  frequently  chosen  to  have  finite  support).  Careful 
attention  is  given  in  the  literature  to  the  rates  of  decay  of  ip]  in  particular,  the  standard  polynomial- 
reproduction  /  quasi-interpolation  argument  requires  ip  to  decay  at  oo  at  a  faster  rate  as  the 
attempted  approximation  order  increases. 
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For  the  purposes  of  this  section,  we  require  the  localization  process  to  satisfy  the  following 
three  conditions: 

Localization  Conditions  2.1.2.  Given  (f>,  and  ip  as  above,  we  assume  throughout  this  section 
that: 

(aj  For  some  >  d, 

(2.1.3)  l/i(a)|  =  Odal"’”'*),  as  joj oo. 

(b)  For  some  >  d, 

(2.1.4)  |V'(^)|  =  <5(!2:|“"**),  as  |x|  -♦  oo. 

(c)  The  sum  in  (2.1.1)  converges  uniformly  on  compact  sets. 

Several  comments  concerning  the  above  assumptions  are  in  order.  First,  the  above  three 
conditions  should  be  regarded  as  very  mild,  and  are  satisfied  by  all  specific  localization  processes 
that  we  are  aware  of.  Second,  in  many  cases  condition  (c)  imposes  decay  rates  on  /i  higher  than 
those  required  in  (a),  making  the  latter  condition  redundant.  Finally,  condition  (c)  is  not  essential, 
and  other,  sometimes  weaker,  notions  of  convergence  can  be  assumed. 

Thanks  to  the  decay  properties  of  ip  at  oo,  it  is  easier  to  write  explicit  approximation  maps  in 
terms  of  its  shifts  than  in  terms  of  the  shifts  of  <p.  Further,  almost  all  analyses  of  approximation 
orders  require  the  basis  function  to  decay  at  oo.  In  general,  a  linear  approximation  scheme  L  which 
explicitly  employs  the  shifts  of  ip  has  the  form 

(2.1.5)  L.f^Y.  ^(•-«)(A/)(a), 

with  /  ►-*  {(A/)(a)}Qg2d  some  linear  assignment.  At  a  minimum,  one  assumes  that  A  is  a  bounded 
map  from  Loo(IR.‘^)  H  C(1R‘^)  (equipped  with  the  uniform  norm)  to  foo(2Z'^).  In  such  a  case,  the 
decay  assumption  imposed  on  ip  in  Localization  Conditions  2.1.2  guarantees  the  sum  in  (2.1.5)  to 
converge  uniformly  on  compact  sets.  It  is  very  convenient  to  assume  further  that  A  commutes  with 
integer  shifts,  i.e.,  that  A/(a)  =  A(/(.  +  a))(0).  For  our  purposes,  it  is  important  to  assume  slightly 
more:  we  need  A  to  commute  with  differentiation  or,  equivalently,  to  be  a  convolution  operator: 

Uniform  Scheme  Condition  2.1.6.  We  assume  that  the  approximation  map  L  is  of  the  form 

~  E  <!.(•- a)  (A/)(a), 

with  A  a  convolution  operator  which  is  well-deSned  and  continuous  from  the  subspace  LooOR’^)  H 
G(IR‘*)  of  LooOR'^)  into  itself. 

As  a  matter  of  fact  (cf.  §3),  in  the  examples  that  we  consider  here,  we  have  either  A  =  identity, 
or  A  =  A*  for  some  band-limited  Li-function  A. 

2.2.  The  map  La  and  its  corresponding  approximation  orders 

With  L  as  above,  we  define  our  approximation  operator  La  that  uses  {(p{-  -  C)}€6H  as  follows. 
First,  recall  from  the  introduction  that  for  each  a  €  ZZ**  we  have  available  a  “pseudo-shift”  <pa  of 
the  form 

(2.2.1)  <Pa  =  '£A{a,0<f>{--0- 

<es 
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At  the  present  stage,  we  do  not  discuss  the  exact  type  of  the  convergence  in  the  definition  of 
<l>a.  Moreover,  our  error  analysis  would  not  make  any  use  of  the  fact  that  ipa  vs  a.  combination  of 
scattered  shifts  of  <f>,  but  only  of  the  fact  that  <j>a  is  “close”  in  some  sense  to  0(-  -  a).  Still,  of 
course,  expressing  <^q  as  a  combination  of  the  appropriate  scattered  shifts  of  <f>  is  pertinent  to  the 
basic  idea  of  approximating  by  scattered  shifts  of  <p. 

The  new  pseudo-shifts  {<Pa}a&‘S,*  localized,  and  by  the  same  localization  sequence  /x 

that  is  used  in  the  uniform  case.  Precisely,  we  define 

(2.2.2)  :=  ^ 

and  substitute  these  “localized  pseudo-shifts”  for  the  shifts  of  ip  in  the  definition  of  L: 

(2.2.3)  Definition  of  La-  With  <p,  ip,  L  and  {ipa}  as  above,  we  define  the  scattered-center  variant 
La  of  L  by 

(2.2.4)  La-I  ^  52  A/(o), 

ae2Z‘‘ 

with  A  the  same  convolution  operator  that  is  used  to  define  L. 

Prior  to  any  analysis  of  the  scheme  La,  we  need  to  show  that  {V'ala  and  La  are  well-defined, 
i.e.,  that  the  infinite  series  used  in  their  definitions  converges  in  some  topology.  This  is  done  with 
the  aid  of  the  following  additional  assumption,  which  is  at  the  heart  of  our  approach: 

Central  Condition  2.2.5.  For  some  mA  >  d, 

Ms  -  .>(•  -  /3))(i)|  <  c(l  +  |i  -  p  6  2Z-. 

with  c  independent  of  x  and  0. 

We  stress  that  none  of  the  results  of  this  section  require  the  representation  (2.2.1),  but  only 
the  above  Central  Condition. 

Invoking  this  new  condition,  we  settle  in  the  next  lemma  the  question  of  the  well-definedness 
of  La- 

Lemma  2.2.6.  Assume  that  the  Localization  Conditions  2.1.2  and  the  Central  Condition  2.2.5 
hold.  Then,  for  a  G  TL^, 

52  ~  /3))(®)1  ^  c'(l  -I-  |i  -  qD""*', 

with  m!  :=  min{m^,mA},  and  d  some  a-independent  constant.  Also 

(2.2.7)  It/-a(x)|<c"(l  +  lx-al)-"*", 

with  m"  :=  xnin{m,j,,mft,mA}  and  c"  some  a-independent  constant. 

Proof.  We  consider  the  difference 

l^a(x)  -  ip{x  -  a)|  <  W  -  a)(</'/?(x)  -  <p{x  -  /3))|. 
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Because  of  Condition  2.2.5,  the  above  sum  is  majorized  by 

(2.2.8)  c  WA- 0)1(1 +  =c  53  W«l(l +  !(!(  + o') -;J|)-”•^ 

&^TL* 

with  l^loo  <  1/2,  and  a'  €  7L^.  By  (b)  of  Localization  Conditions  2.1.2,  and  for  any  fixed  y,  this 
latter  sum  is  the  evaluation  at  a'  of  the  discrete  convolution  of  sequences  which  decay  at  oo  like 
the  m^th  and  ru/ith  inverse  power.  Since  y  is  restricted  to  a  compact  set,  there  exist  constants 
such  th'^t 

ci(i  +  \y-  tI)"'”'*  <  (1  -t-  W)"'"'*  <  C2(l  +  12/  -  71)”'"'*. 

which  allows  us  to  bound  the  sum  in  (2.2.8)  by 

const(l  +  |a'|)""*  <  const'(l  +  I2/ +  a'l)""*  =  const'(l  +  |a;  -  , 

for  m'  :=  min{m^,m>i},  and  for  some  a-independent  constants  (cf.  e.g.,  the  argument  in  the 
proof  of  Proposition  6  in  [BuD]).  Combining  this  with  (b)  of  Conditions  2.1.2  we  obtain  the  decay 
assertion  onxpa-  ^ 

We  are  now  ready  to  state  and  prove  the  main  theorem  of  this  section: 

Theorem  2.2.9.  Let  cf),  n,  ipi  and  L  be  as  above,  and  assume  that  the  Localization  Conditions 
2.1.2,  and  the  Uniform  Scheme  Condition  2.1.6  are  satisfied.  Let  ((t>a)aez<‘  ^  collection  of 

functions  that  satisfy  the  Central  Condition  2.2.5,  and  let  La  be  their  associated  map  defined  as 
in  (2.2.4).  Assume  further  that  for  some  positive  integer  k  and  0  <  e  <  1  the  following  hold: 

(a)  )/i(a))  =  as  lal  -♦  00. 

(b)  The  functional  fi  defined  by 

Mp:= 

ae.Tl* 

annihilates  Ilfc. 

Then,  for  every  f  €  n 

11(1  -  1a)(ta/)|U  <  constd/U.,  +  I/U.I+,)  I  Jt+l]  J  ^  J; 

Furthermore,  ife  =  1,  and  (a)  above  is  replaced  by  the  stronger  condition  Eaez"  <  00. 

then,  for  f  as  above, 

||(L  -  La)(o-/./)I|oo  <  const  I/loo, fc+i- 

Proof.  We  first  prove  the  last  (and  simplest)  case  stated  in  the  theorem,  i.e.,  when 
Eaez'*  <  OO-  will  then  snow  how  to  modify  the  argument  of  that  case  to  ob¬ 

tain  the  results  stated  with  respect  to  the  other  cases. 

Given  /  €  W^'*'*(1R‘^)  n  ®  €  IR*^,  we  first  estimate 

M/(a:  +  -))=  M-a)/(a:  +  Of)- 

ae2i<‘ 

With  Txf  the  Taylor  expansion  of  degree  k  of  f  about  x,  the  fact  that  p  annihilates  Ilfc  implies 
that  p{Txf)  =  0.  Since  |/(x -Ha)  -  (Txf){x  -H a)|  <  const|/|oo,*+i  |a|*''*'^  we  obtain  the  estimate 

(2.2.10)  \p{f{x  +  •))!  <  const(/|oo,fe+i  =  const'|/|oo,)fc+i. 

aeZZ'* 
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Now,  we  can  write  (L  -  LA){(Thf)  as  follows: 


1(1  -  LA)i<rhf)(x)\  =1  ^  (V'(-  ~a)~  V’a)(^)  A(cr/,/)(a)| 

=1  52  -fi)-  4>0{x))fJ.{P  -  a)  A(tT/,/)(a)l. 

a€2Z^ 

By  Lemma  2.2.6  and  Condition  2.1.6,  the  above  double  sum  converges  absolutely,  and  summation 
by  parts  yields  that 

(2.2.11)  \{L  -  La){ctMx)\  =  1  52  M^))  E  -  ‘>)A(^Ta/)(a)|. 

/Sea''  ae7Z‘‘ 

Next,  we  identify  the  sum  y-iP  -  Q)i\{ahS){ot)  in  (2.2.11)  as  /i(A(cr/i/)(/3  +  •)))  and 

invoke  (2.2.10)  to  estimate  the  latter  expression  as  follows: 

(2.2.12)  l/x(A(cTh/)(/3  +  •))!  <  constjA  *  {(Thf)\oo,k+i  <  const  h'''*'^||A|||/|oo,fc+i, 

where,  in  the  second  inequality,  the  fact  that  convolution  commutes  with  differentiation  and  that 
dilation  is  an  isometry  on  Loo(ni‘^)  are  used.  Substituting  (2.2.12)  into  (2.2.11),  we  arrive  at  the 
bound 

|(L  -  L^)(cT/i/)(z)i  <  const  h'‘'*'^\f\oo,k+i  E 

I}e2z<‘ 

which  completes  the  proof  of  the  present  case,  since  Condition  2.2.5  implies  that  the  last  sum  is 
majorized  by  a  (periodic)  bounded  function. 

In  case  (a)  and  (b)  hold,  but  the  stronger  assumption  53o€Z''  <  oo  is  not  valid, 

we  modify  the  argument  that  leads  to  (2.2.12)  as  follows.  We  partition  IR'^  to  a  ball  Bh  of  radius 
1/h  centered  at  the  origin  and  to  the  complement  of  that  ball.  If  a  6  B/,,  we  estimate  the 
difference  f{x  +  q)  -  {Txf){x  +  a)  as  before,  i.e.,  obtain  that 

|/(x  +  Of)  -  (Tx/)(i  +  a)|  <  const|/|oo,*+i|al'''^^ 

However,  for  a  €  B^,  we  treat  the  zero  f{x  +  •)  “  {TxJ){x  +  •)  has  at  the  origin  as  of  order  k 
(though,  it  is  of  order  k  +  1).  This  leads  to  an  estimate  of  the  form 

\f{x  +  a)  -  (Tx/)(z  +  a)|  <  const|/|oo,fclal''. 

Summing  the  first  estimate  over  a  6  B/i  and  the  second  one  over  a  G  and  invoking  the  decay 
assumption  (a)  on  y,  we  get  that 

IM/(x  +  .))|<const(|/|oc.it+i  E  +  E  H'""')- 

o€Bi.  oeB' 

According  to  Lemma  4.2  of  [DJLR],  the  sum  in  the  second  summand  above  behaves  like  0{h^), 
while  the  first  one  is  0{h^~^)  and  0(|  log  h|)  for  0  <  e  <  1  and  e  =  1,  respectively.  We  thus  obtain 
an  estimate  analogous  to  (2.2.10),  which,  say  for  e  <  1  reads  as 

lM(/(a:  +  •))!  <  const(|/|oo.fc+i/i'“‘  +  |/|oo,fcA'). 
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Following  the  argument  before  (2.2.12),  we  obtain  that 

l/*(A(ff/i/)(/3  +  •))!  <  const  h''+'||All(l/loo,fc+i  +  l/loo,fc). 

The  proof  is  now  obtained  as  in  the  previous  case,  with  (2.2.12)  being  replaced  by  the  above  bound. 
The  case  e  =  1  is  similar.  4 

At  a  first  look,  it  is  not  clear  that  the  theorem  allows  us  to  draw  conclusions  on  the  approx¬ 
imation  orders  provided  by  the  map  La,  and  because  of  two  reasons:  first,  the  theorem  assumes 
that  fi  annihilates  Ilfc,  and  one  may  suspect  this  further  condition  on  /x  to  compete  with  other 
requirements,  or,  at  least,  to  exclude  various  localization  processes  which  do  not  prepare  for  such 
a  condition.  Second,  the  powers  and  that  appear  in  the  conclusion  of  the  theorem  seem 
unrelated  to  the  possible  approximation  order  provided  by  the  map  L.  The  next  two  results  are 
meant  to  clarify  these  two  important  points.  The  first  shows  that  the  annihilation  assumption  on  fi 
is  usually  a  by-product  of  the  Localization  Conditions  2.1.2.  The  second  shows  that  in  most  cases 
of  interest  the  parameter  k  +  £  that  appears  in  Theorem  2.2.9  is  no  smaller  than  the  approximation 
order  provided  by  L. 

Lemma  2.2.13.  Let  k  be  a  positive  integer,  and  assume  that  (p  is  continuous  on  1R‘*\0  and  has  a 
singularity  of  order  >  k  at  the  origin.  Assume  further  that  the  Localization  Conditions  2.1.2  hold, 
and  that  the  linear  functional 

^  M-a)p(a) 

is  well-defined  on  Uk  (i-e.,  the  above  sum  converges  absolutely  for  every  p  G  n,fc).  Then  p  annihilates 

Uk. 

Proof.  Condition  (c)  of  Localization  Conditions  2.1.2  allows  us  to  compute  the  Fourier 
transform  of  ^  in  (2.1.1)  term  by  term.  Defining 

(2.2.14)  P-= 

Qez** 

we  thus  obtain  that,  pointwise  on  1R‘^\0,  i)  =  4>p.  Due  to  (b)  of  Localization  Conditions  2.1.2,  V'  is 
continuous  everywhere,  and  in  particular  at  the  origin,  and  hence  p  necessarily  has  a  zero  of  order 
>  k  at  the  origin.  Since  p  is  well-defined  on  I!*,  p  is  fc-times  differentiable  at  the  origin  (and,  as 
a  matter  of  fact,  everywhere),  and  can  be  differentiated  term  by  term.  Consequently,  given  any 
p  G  Ilfc,  we  obtain 

0  =p(-zD)p(0)  =  p{a)p{-a)  =  pp, 

a€Z'' 

as  claimed.  A 

Lemma  2.2.15.  Assume  that 

(a)  The  Localization  Conditions  2.1.2  hold; 

(b)  <l>  is  continuous  on  IR  \0,  has  a  singularity  of  order  k'  at  the  origin,  and  satisfies,  for  some 

2ir7L\Q,  $(J)  ^  0; 

(c)  ^(0)  #  0. 

Then  the  approximation  order  provided  by  L  to  W^~(IR‘^)  is  <  k'. 
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Proof.  The  proof  follows  from  the  general  analysis  of  approximation  orders  of  [BR],  and 
the  claim  is  essentially  proved  in  Theorem  4.2  there,  hence  is  only  sketched  here.  First,  Theorem 
2.17  of  [BR]  implies  that  if,  for  some  &  €  IR*^,  L  provides  approximation  order  n  to  the  exponential 
function  e#,  then,  for  any  j  6  27rZZ‘*\0,  and  in  particular  for  the  one  specified  in  (b)  above, 

m  +  he)\  =  o{h^). 

Second,  since  1/V’  is  bounded  in  a  neighborhood  of  zero,  the  last  condition  implies  that 

fO'  +  he)/^(he)  =  (9(/i"). 

Since  tp  =  'ii<p,  and  j  is  a  period  for  p,  we  obtain  that 

^{j  +  he)/Hh9)  = 

Since  (f>{j)  ^  0,  we  finally  conclude  that,  in  case  L  provides  approximation  order  n  for  all  exponen¬ 
tials  €0,  6  €  IR'^, 

i/Pihd)  =  Oik'll  \/e  e  iR'', 

and  since,  by  assumption,  p  has  a  singularity  of  order  k'  at  the  origin,  we  conclude  that  n  <  k'.^ 

The  following  theorem,  which  is  a  direct  consequence  of  Theorem  2.2.9  and  the  last  two  lemmas, 
summarizes  the  results  obtained  for  functions  ^  whose  Fourier  transform  has  a  singularity  at  the 
origin. 

Theorem  2.2.16.  Assume  that: 

(a)  ^  is  continuous  on  1R‘^\0,  and  has  a  singularity  at  the  origin  of  exact  order  k'. 

(b)  The  Localization  Conditions  2.1.2  and  the  Uniform  Scheme  Condition  2.1.6  hold. 

(c)  For  some  integer  k  <  k',  and  e  <1,  n  is  either  well-defined  on  (e  =  1)  or  satisfies 

|/i(Q)|  =  as  |q|  — >  oo, 

(e  <  1). 

(d)  The  Central  Condition  2.2.5  holds. 

Then: 

(i)  La  provides  the  same  approximation  order  as  L  provides  to  functions  in  iy^‘^^(IR‘*)nC*'‘*‘^(]R‘^), 
unless  this  latter  order  exceeds  k  -i-e. 

(ii)  Ifk  +  e  =  k',  La  necessarily  provides  the  same  approximation  order  as  provided  by  L  (to  the 
space  in  (i))  in  case  (p  is  non-zero  at  some  j  €  2irZ  \0,  and  V’(O)  ^  0. 

Remark  2.2.17.  Condition  (c)  in  the  above  theorem  distinguishes  between  an  integral  approxi¬ 
mation  order  and  fractional  approximation  order,  and  is  more  restrictive  in  the  former  case.  Had 
we  merely  assumed  that,  for  s  =  1,  |/x(o)|  =  Theorem  2.2.9  would  have  yielded  an 

approximation  rate  of  logftj),  rather  than  the  obtained  above.  Though  there 

are  situations  (see  §3)  where  only  the  slower  decay  rate  is  known  about  p,  the  approximation  rate 
L  is  known  to  provide  (in  all  such  cases  that  we  are  aware  of)  is  log  A|),  so  that  La  still 

inherits  the  approximation  power  of  L. 
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Remark  2.2.18.  Part  (ii)  of  the  last  theorem  deals  with  “the  worst  case”  showing  that,  for  some 
functions  /,  the  approximation  order  provided  by  La  to  /  does  not  lag  behind  the  one  provided 
by  L.  It  does  not  exclude  the  possibility  that,  for  other  functions  /  e  H^4'*'^(1R‘^)  D  L 

provides  approximation  order  n>  k'  to  /,  and  for  such  /,  La  may  fail  to  provide  approximation 
order  n  for  /,  since  it  is  guaranteed  to  provide  only  order  k'.  However,  we  expect  such  a  situation 
to  be  truly  exceptional,  and  more  precisely  we  conjecture  that,  under  the  assumptions  of  Theorem 
2.2.16,  no  function  /  €  n  other  than  the  constants  can  be  approximated  to  a 

rate  better  than  k'.  In  this  regard,  we  mention  that  (R)  proves  that  in  the  1,2-analog  of  the  present 
situation  no  non-zero  function  in  the  potential  space  IV2  can  be  approximated  to  an  order 

>ib'. 

2.3.  A  complementary  error  analysis  of  the  approximation  order  of  La 

Theorem  2.2.16  asserts  that,  under  various  conditions,  the  approximation  order  provided  by 
La  is  no  smaller  than  the  one  provided  by  L.  Most  of  the  assumed  conditions  are  acceptable  in 
the  sense  that  they  hold  in  all  the  examples  that  are  analysed  in  the  next  section.  However,  one  of 
the  conditions  assumed  in  Theorem  2.2.16  may  appear  to  be  too  restrictive,  and  excludes  several 
examples  of  interest;  we  assume  that  /x  decays  fast  enough  to  be  absolutely  summable  against 
polynomials  in  H*;.  This  condition  can  be  relaxed  by  using  Fourier  analysis  alternatives,  as  we 
discuss  in  the  present  subsection. 

In  the  next  theorem  we  use  the  notation 

to  denote  the  space  of  all  functions  /  that  (i):  their  Fourier  transform  /  is  a  Radon  measure,  and 
(ii):  the  total  mass  ||(1  -f  |  •  l)*/l|i  of  (1  -1-  ]  •  |)*/  is  finite.  The  norm  |1(1  -t- 1  •  D^/K^oincides  with 
the  Z/i(lR‘^)-norm  of  (1  -I- 1  •  ()*'/,  whenever  /  is  a  function.  It  induces  a  norm  on  which 

we  denote  by  ||  •  ||^  *,  i.e., 

- 11(1  + 1  • 

Similarly,  we  use 

l/L.* 

to  denote  the  total  mass  |||  •  l*/|ji  of  |  •  j*'/.  Note  that,  for  integer  A:,  VF^(1R‘^)  is  continuously 
embedded  into  1F’^(1R‘^)  fl  C*(IR“). 

Theorem  2.3.1.  Assume  that  Centra]  Condition  2.2.5,  Localization  Conditions  2.1^,  and  Uniform 
Scheme  Condition  2.1.6  bold  and  that  the  convolution  function  A  of  2.1.6  satisfies  A  6  I,oo(IR‘*)- 

(a)  If,  for  some  positive  k,  |  •  is  bounded,  then,  for  every  f  £  W^(1R'^), 

||(L  -  LA)i<Thf)\\oo  <  const 

In  particular, 

(b)  La  provides  to  functions  in  1F^(1R‘^)  the  same  approximation  order  provided  by  L  in  case  <t>  is 
continuous  on  IR‘^\0,  has  a  singularity  of  exact  order  k'  at  the  origin,  does  not  vanish  at  some 
j  €  2^22*^,  and  ip(0)  ^  0. 

Proof.  We  first  claim  that: 
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Lemma.  Under  the  conditions  of  part  (a)  of  Theorem  2.3.1, 


\fig\  <  constl^r^  jk,  g  £ 


Proof  of  the  Lemma.  We  first  assume  that  g  £  is  compactly  supported,  and  denote 

^1  •“  Then,  with  5|  the  Fourier  series  of  g^, 

(2.3.2)  pg=  p{-a)g{a)  =  ((i*g\)(0)  =  (27i)~‘'  f  il{t)g^it) dt. 

aez** 

Since  g  £  W^(IR‘^)  and  is  compactly  supported,  g  £  Li(lR‘*),  and  therefore  we  may  invoke  Poisson’s 
summation  formula  to  obtain  that 

^  ?(■  +  a)  =  (2<r)-'  Yi  9Wy-e  =  (2ir)-'j|. 

a62irZ‘‘ 


Combining  the  above  with  (2.3.2),  we  arrive  at 


(2.3.3)  P9=  f  Ki)  ^  git  +  a)dt=  ^  f  Ji{t)g{t  +  a)  dt  =  f  pg. 

Q€22rZZ‘‘ 


a62r2‘‘ ' 


(The  summation  by  parts  in  the  second  equality  can  be  justified  by  dominated  convergence  argu¬ 
ments,  since  p  is  bounded  and  the  sum  there  is  Li-convergent.) 

Next,  we  extend  (2.3.3)  to  a  general  g  £  kF'^(lR‘^).  For  that,  we  take  {T]h)h  to  be  a  band-limited 
rapidly  decaying  approximate  identity,  with  rji  equals  1  around  the  origin,  and  define  gh  :=  gv^- 
By  (2.3.3), 

m  =  /  MS/i  =  (27r)“‘'  f  p{g  *  r]h)  =  (27r)"  V  *  (5  ♦  p{—)){0). 

Jsi<‘  Jti<^ 

Clearly,  pg^  —*  (2Tr)~*^ pg  as  h  0.  Also,  since  g*p{ — )  is  continuous  and  bounded,  its  convolution 

with  the  approximate  identity  converges  pointwise  to  itself,  and  hence 

Vh*{9*fi{—))i0)^{g*p{—))i0)=f  pg,  as  h -*  0, 

Jr* 

and  the  desired  extension  of  (2.3.3)  easily  follows. 

Since,  by  assumption,  |^|  <  const]  >  I*",  we  can  deduce  now  that 

ImsI  =  I  5^  /*(-a)<7(«)l  <  const  f  |  •  1*  |^  =  constl^l'oo.fc. 

4k 

We  now  continue  with  the  proof  of  Theorem  2.3.1.  Similarly  to  the  proof  of  Theorem  2.2.9, 
we  estimate  p{A.{ahf){P  +  •))•  Invoking  the  lemma  with  g  :=  A.{ahf){l3  +  •),  we  obtain  from  the 
translation  invariance  of  the  semi-norm  |  •  |(,j,  ^  that 

lM(A(<Tfc/)(/3  +  •))!  <  const|A  *  (<Th/)iJ,o.fc- 
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Since  A  is  bounded,  A*  can  be  easily  seen  to  be  continuous  on  W'^(]R‘‘),  and  hence,  by  elementary 
properties  of  Fourier  transform 

+  •))!  <  const|<T/i/l'^  t  =  const^i^l/l^o,*- 

The  rest  of  the  proof  of  (a)  is  identical  with  its  counterpart  in  Theorem  2.2.9. 

In  the  proof  of  (b),  we  employ  Lemma  2.2.15.  Note  that  the  lemma  is  stated  with  r^pect 
to  IV^(1R‘*)  n  C*'(]R‘*),  but  is  valid  also  with  respect  to  the  smaller  space  IV^(l^),  since  the 
exponential  functions  6$,  6  €  that  are  used  in  the  proof  there  are  contained  in  too. 

Under  the  further  eissumption  in  (b),  we  can  use  an  argument  as  in  the  first  part  of  the  proof  of 
Lemma  2.2.13  to  conclude  that  /I  has  a  zero  of  order  k'  at  the  origin,  and  therefore  the  conditions, 
hence  the  consequences  of  part  (a)  here  are  valid,  for  the  choice  k  :=  k'.  Combining  part  (a)  with 
Lemma  2.2.15,  we  obtain  (b).  ♦ 


2.4.  Approximation  to  unbounded  functions 

The  kind  of  approximation  which  is  considered  in  the  previous  section  is  stationary.  That 
means,  by  definition,  that  the  space  from  which  our  approximant  to  is  selected  does  not  change 
with  h.  On  uniform  grids,  stationary  approximation  have  the  following  important  advantage:  in 
many  cases  the  approximation  map  L  can  be  extended  in  a  natural  way  to  a  map  L  that  acts  from 

(2.4.1)  Nk  :=  {/  :  IR''  ^  C  :  1|/|U*  l!(l  +  1  •  |)~Vllt«(R-)  <  oo) 

into  S^d{<f)),  and  still  provides  to  A:-time  continuously  differentiable  functions  in 

<(1R'‘)  :=  {/  €  AT*  ;  <  cx.) 

the  same  approximation  order  as  provided  by  L  to  lV^(]R'')nC‘(lR'*).  If  L  has  the  form  assumed 
in  Uniform  Scheme  Condition  2.1.6,  i.e., 

(2.4.2)  53  t^(-  -  o)(A/)(a), 

the  “natural  extension”  L  still  satisfies  the  same  rule,  and  one  then  needs  to  ensure  that  A*  is 
well-defined  and  continuous  on  Nk,  and  that  the  sum  in  (2.4.2)  is  meaningful  (i.e.,  converges  in 
soir**  reasonable  sense)  even  for  f  £  Nk- 

In  case  L  is  well-defined  and  is  known  then  to  provide  approximation  order  k  not  only  to 
ifc-time  differentiable  functions  in  IU^(Bl'*),  but  also  in  the  larger  space  it  is  desirable  to 

extend  La  appropriately.  However,  for  such  an  extension,  we  need  to  assume  slightly  more  than  the 
assumptions  made  in  Theorem  2.2.9.  The  modifications  are  mainly  concerned  with  the  imposition 
of  faster  decay  rates  on  (i  and  ^  in  Localization  Conditions  2.1.2,  and  on  (j)c,  —  <!>{•  —  oc)  in  Central 
Condition  2.2.5.  It  is  worth  mentioning  that,  at  least  as  far  as  the  decay  rates  of  ^  are  concerned, 
the  modification  seems  to  be  right  in  place:  for  some  /  €  u;^(lR'*),  we  expect  the  coefficients 
{A(/)(a)}<,g2'‘  grow  1*^0  !“!*'>  hence  need  the  faster  decay  rates  for  the  convergence  of  the  series 

(2.4.3)  1/  :=  53  "  “)(A/)(«)- 
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Theorem  2.4.4.  Assume  that  LocaUzation  Conditions  2.1.2,  as  well  as  Central  Condition  2.2.5 
bold  for 

(2.4.5)  m^,  m^,  >  d  +  r, 

for  some  positive  r.  Assume  further  that  A  of  (2.1.5)  extends  continuously  to  a  convolution  map 
from  Nr  into  itself  Then,  the  series  (2.2.4)  that  debnes  La/  converges  absolutely  and  uniformly 
on  compact  sets,  for  every  f  €  Nr.  Further,  assume  that  conditions  (a)  and  (b)  in  Theorem  2.2.9 
bold,  with  k  there  >  r.  Then,  the  assertions  made  in  Theorem  2.2.9  about  {L  -  LA){<Thf)  hold  for 
functions  in  n  O  H  Nr. 

Remark.  The  standard  choice  for  r  above  hr  =  k.  The  fact  that  we  allow  here  r  <  k  as  well,  is 
aimed  at  deriving  partial  improvements  of  Theorem  2.2.9  for  the  case  when  the  decay  rates 
are  only  partially  improved. 

Proof.  Since  we  are  assuming  that  >  d  +  r,  we  conclude  from  Lemma  2.2.6 

that 

(2.4.6)  |V'a(a:)l  <  const(l  +  |x  -  al)"*” 


with  m  >  d  +  r,  and  since,  further.  A/  €  Nr,  the  map  La  is  well-defined  on  Nr- 

We  may  then  follow  the  proof  of  Theorem  2.2.9  verbatim,  with  straightforward  modifications 
of  some  of  the  justifications  of  the  various  steps.  For  example,  the  fact  that  dilation  is  an  isometry 
on  Ioo(IIl‘*)  is  still  effective  in  (2.2.12),  since  the  estimate  of  |A  *  <Thf\oo,k  involves  the  convolution 
of  A  with  bounded  functions  only  (despite  of  the  fact  that  /  itself  is  not  necessarily  bounded).  As 
another  instance,  the  summation  by  parts,  which  is  applied  to  the  double  sum  in  the  (2.2.11),  is 
also  justified  by  the  modified  decay  assumptions  on  |^(-  -  0)  ~  <p0\  combined  with  the  fact  that 
A(<r„/)  €  Nr.  ♦ 

Theorem  2.4.4  leads  to  an  improvement  of  Theorem  2.2.16  by  extending  the  range  of  approx- 
imants  from  ^^+‘(5^'')  H  to  p,  ti;^(lR'*)  n  C^+‘(IR'')  n  Nr.  We  omit  these 

straightforward  details. 


2.5.  Polynomial  reproduction 

With  ^  as  in  Localization  Conditions  2.1.2,  and  L  as  in  (2.1.5),  we  say  that  L  reproduces  Ilr 
(r  >  0),  if  is  the  identity.  Embedded  in  this  is  the  requirement  that  of  Conditions  2.1.2 
would  be  no  smaller  than  a  4  r.  The  polynomial  reproduction  property  is  an  essential  ingredient  in 
the  error  analysis  of  the  so-caiied  quasi-interpolation  schemes,  and  in  turn,  makes  the  approximation 
orders  that  can  be  established  by  such  schemes  restricted  by  the  decay  rate  of  ip. 

In  our  context  here,  we  do  not  need  to  know  the  argument  used  for  the  derivation  of  the 
approximation  order  of  L,  and  more  significantly,  do  not  need  a  corresponding  polynomial  repro¬ 
duction  property  from  La.  This  is  very  important  because  of  the  following:  if  L  is  proved  to 
provide  approximation  order  k  via  a  polynomial  reproduction  argument,  then  we  know  that  the 
localization  ip  employed  by  L  satisfies  >  d  -i-  k  —  1.  In  order  to  show  that  La  also  reproduces 
polynomials,  we  need  to  know  that  La  is  well-defined  on  Ilj^-i,  which,  if  Lemma  2.2.6  is  used  for 
that,  will  require  further  assumptions,  primarily  that  of  Condition  2.2.5  would  satisfy,  at  a 
minimum, 

>0  +  ^-1. 
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This,  by  all  means,  is  not  a  minor  requirement,  since  it  makes  the  construction  of  the  matrix  A 
more  involved  (see  §2.7).  Fortunately,  even  though  the  approximation  orders  of  L  might  have  been 
derived  via  polynomial  reproduction,  our  conversion  results  (Theorem  2.2.9  and  Theorem  2.2.16)  do 
not  require  La  to  reproduce  polynomials,  and  hence  we  are  not  forced  to  look  for  A  that  improves 
upon  the  basic  condition  >  d. 

However,  it  is  worth  mentioning  that  the  polynomial  reproduction  property  of  L  is  inherited 
by  La  if  the  various  decay  rates  m^,m^,myi  are  known  to  be  large  enough: 

Corollary  2.5.1.  Under  the  notations  and  assumptions  of  Theorem  2.4.4,  if  L  reproduces  Hr 
(r  <  k)  so  does  La- 

Proof.  Let  p  €  Hr.  Had  we  known  that  degp  <  k,  the  conclusion  would  have  been 
immediate,  since,  by  Theorem  2.4.4,  ||(£  -  La)p\\  <  const(|p|oo,fc  +  bloo,fe+i)  =  0.  The  argument, 
however,  for  the  case  degp  =  k,  is  not  much  more  involved:  since  we  assume  that  p  annihilates 
Hfc,  hence  Hr  (that  assumption  is  made  in  Theorem  2.2.9,  and  is  adopted  in  Theorem  2.4.4,  hence 
here),  we  may  derive  from  the  proof  of  Theorem  2.2.9  (cf.  (2.2.11)  and  its  following  line)  that 

(2.5.2)  (L  -  La)p  =  +  0))- 

Since  A  is  convolution.  Hr  is  an  invariant  subspace  of  it,  and  therefore  Ap  €  Hr-  Since  p  annihilates 
Hr,  (2.5.2)  implies  that  Lap  =  Lp  =  p,  the  second  equality  by  the  polynomial  reproduction 
assumption  made  with  respect  to  L.  ^ 

The  crucial  assumptions  in  the  above  corollary  are  (a)  and  (b)  of  Theorem  2.2.9,  stating  that 
p  is  well-defined  on  Hr  and  annihilates  it.  It  is  possible  to  replace  thesejwo  assumptions  by  the 
condition  >  d  +  r,  if  we  assume,  as  we  do  in  the  present  paper,  that  <(>  is  continuous  on  1R‘*\0, 
and  that  ^  does  not  vanish  identically  on  27rZZ‘*\0.  Indeed,  it  is  well-known  then,  that  th^Hr- 
reproduction  property  of  L  implies  that  |  •  (f>  is  singular  at  the  origin,  which  means  that,  if  <f>  has 
a  singularity  of  some  order  at  the  origin,  that  order  is  >  r.  Consequently,  Lemma  2.2.13  implies 
that  p  annihilates  Hr. 

The  argument  used  in  the  proof  of  Corollary  2.5.1  supports  also  the  following  stronger  claim. 

Corollary  2.5.3.  Let  P  be  a  translation-invariant  space  of  polynomials,  which  is  reproduced  by 
L  and  is  annihilated  by  p.  If,  further,  the  double  sum 

a.PeZ"' 

(cf.  (2.2.11))  converges  absolutely,  for  every  p  G  P,  then  La  reproduces  P,  as  well. 

2.6.  Remarks  on  the  location  of  the  centers 

In  actual  approximation  schemes,  we  do  not  dilate  the  approximand  /.  Instead,  we  dilate  the 
space  Ss((f))  from  which  the  approximants  are  selected.  Thus,  the  “true”  error  to  be  considered 
should  have  been 

/  -  <T\lhL{(Thf)- 

The  change  to  the  study  of  a/,/  -  L{<Thf)  is  due  to  technical  convenience  and  is  made  available  by 
the  invariance  of  the  max-norm  under  dilation. 
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However,  in  order  to  understand  the  nature  of  our  approximants  as  h  changes,  it  is  instructive 
to  look  closer  at  the  approximant  <Ti/hL{ahf)  and  its  scattered  center  analogue.  First,  we  see  that, 
adopting  the  convolution  assumption  on  A  (cf.  Uniform  Scheme  Conditions  2.1.6), 

<ri/hL{<Thf)  =  ^  (•  -  Q!)(Afc/)(Q), 

where 

Afc  :  /  hr^{ciiHX)  *  f. 

The  scattered  center  version  OxihLA{.Ohf)  is  obtained  by  replacing  each  -  a)  by  eri/hi’a- 

Finally,  the  function  o’l/hV'a  is  obtained  by  an  application  of  the  same  matrix  A  to  the  scattered 
shifts  of  only  that  the  center  set  E  is  replaced  by  the  scaled  set  /lE.  Thus,  in  summary,  the 

scattered  center  approximant  aijt^LAiohf)  makes  use  of  the  AE-translates  of  criih<t>. 

The  fact  that  at  the  A-level  the  dilated  center  set  AE  is  employed  should  not  be  regarded  as  an 
essential  ingredient  of  our  approach  and  arguments.  As  a  matter  of  fact,  one  might  employ  at  the 
A-level  a  center  set  E/,  which  may  resemble  no  relation  to  AE.  This,  in  turn,  forces  the  search  of 
different  matrices  Ah,  for  different  values  of  A.  Applying  Ah  to  one  obtains  a  set  of  functions 

and  the  Central  Condition  2.2.5  should  then  be  modified  to 

-  Of)!  <  const(l  +  \x  -  al)""*'*,  a  e  7Z‘^, 

with  const  independent  not  only  of  a  and  x  but  also  of  A. 

Another  remark  is  concerned  with  the  stationary  nature  of  the  approximation  schemes  dis¬ 
cussed  in  this  section.  That  notion  refers  to  the  fact  that  the  sequence  of  approximants  {L{(Thf)}h 
to  {(Thf}h  employs  translations  of  the  same  function  at  all  A-levels.  There  are  important  sit¬ 
uations  in  spline  theory  when  such  simplification  in  the  construction  is  unacceptable  (primarily 
exponential  box  splines,  cf.  [DR]).  But,  in  the  context  of  radial  basis  function  approximation  the 
only  presently  known  case  when  non-stationary  schemes  are  used  is  concerned  with  spectral  approx¬ 
imation  (cf.  [BeLi],  [BuD],  [BR],  [R]),  a  case  that  cannot  be  covered  by  the  approach  developed 
here. 


2.7.  On  the  Central  Condition  2.2.5 

Condition  2.2.5  is  the  major  condition  that  is  assumed  in  both  Theorem  2.2.9  and  Theorem 
2.2.16.  As  a  matter  of  fact,  any  set  {<^a}a  that  satisfies  Condition  2.2.5  can  be  used  to  define 
the  approximation  map  La-  However,  since  our  interest  is  in  approximating  from  the  span  of  the 
scattered  translates  of  <}>,  we  certainly  restrict  our  attention  to  functions  {<^a}  which  are  expressed 
in  the  form  (2.2.1). 

In  the  present  subsection,  we  briefly  discuss  the  following  important  (and  natural)  question: 
what  sufficient  conditions  on  the  coefficient  matrix  A  imply  Condition  2.2.5?  Here  we  show  that, 
under  some  basic  assumptions  that  we  adopt  here  with  respect  to  the  structure  of  the  Fourier 
transform  of  4>,  one  should  aim  at  approximating  (around  the  origin)  the  constant  function  1  by 
exponentials  of  the  form 

Ca  = 

for  each  a  G  7L^.  The  (formal)  derivatives  of  e*  are  of  the  form 
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and  therefore,  given  an  integer  j,  if  the  sequence  $  (1  +  [a  -  4p)i4(a,0  is  in  fi(E:),  then  e* 

is  j-times  boundedly  differentiable,  can  be  differentiated  term  by  term,  and  its  H^^(IR'*)-norm 
satisfies 


Theorem  2.7.1.  Let  A  =  (A(a,0)jz<<xi  coefBcient  matrix  of  (2.2.1),  and  assume  that  each 

of  the  rows  of  A  is  in  Let  e*  be  defined  as  above.  Suppose  the  following  conditions  hold: 

(a)  The  set 

{H  9  ^  ^  (1  +  1^  -  :  a  €  Z*'} 

lies  in  fi(E)  and  is  bounded  there,  for  all  j  <  s,  for  some  non-negative  integer  s. 

(b)  For  each  a  e  S'*,  all  derivatives  of  1  —  e*  of  orders  <  m  vanish  at  the  origin,  for  some 

non-negative  integer  m  <  s.  ^  ^ 

(c)  For  yet  another  positive  integer  r,  (p  is  r-times  differentiable  on  ]R‘*\0,  and  each  D°<p,  |a|i  <  r, 
is  summable  around  oo. 

(d)  Each  D°‘$,  |q|i  <  r,  (calculated  on  IR‘'\0j  has  a  singularity  of  order  k'  +  [ali  at  the  origin,  for 
some  positive  k'  ;=  k'((p),  and  for  the  same  r  of  (c). 

(e)  The  distribution  (p  is  of  order  <  m. 

Then,  Condition  2.2.5  is  satisfied  with  mA=m-  (fc'j  +  d,  with  [fc'j  the  greatest  integer  <  k',  and 
provided  that  this  exceeds  neither  r  nor  s. 

Roughly  speaking,  r  can  be  regarded  as  a  smoothness  parameter  for  (p.  Thus,  in  these  terms, 
the  threshold  for  the  “order  of  matching”  m  (in  (b)  above),  in  order  to  achieve  >  d,  is  the 
order  of  the  singularity  of  ^  at  the  origin.  Above  that  threshold,  a  better  matching  (i.e.,  higher  m) 
would  result  in  a  larger  in  Condition  2.2.5  (which  might  be  important,  e.g.,  if  we  approximate 
unbounded  functions,  cf.  §2.4),  but  overmatching,  i.e.,  matching  which  exceeds  r  +  [*')  -  d,  may 
not  yield  any  improvement  in  ttia- 

Proof.  Condition  2.2.5  is  equivalent  to  the  statement  that  {(1  +  |  •  |)”*^(<^  -  (Pa('  +  oi))  '• 
a  6  Z**}  is  bounded  in  Loo(lR*^)-  Instead,  we  prove  that  the  Fourier  transforms  of  these  functions 
are  bounded  in  L^Hl*^).  Since  =  m  -  {*']  +  d,  it  suffices,  for  that  latter  purpose,  to  prove 
that,  for  each  integer  j  <m  -  (*']  +  d,  the  Fourier  transforms  of  |  •  p(<6  -  d>a{-  +  a)),  a  €  Z**  are 
bounded  in  ii(lR'').  Since  the  Fourier  transform  of  -  <Pa{-  +  Qt)  is  (1  -  6^)0,  we  may  prove  that 
the  functions 

-  e;)^) :  a  e  Z^  W  <m-  [k']  +  d} 

are  bounded  in  Li(lR‘*).  By  applying  Leibnitz’  rule,  we,  finally,  may  replace  this  last  set  of  functions 
by 

(2.7.2)  {£>^(1  -  :  a  €  Z**,  7  <  \P]i  <m-  (fc'l  +  d}. 

Since  m  -  (*']  +  d  <  r,  by  assumption,  condition  (c)  implies  that  D^~'’(p  is  summable  around 
00.  In  addition,  since  m  — (A:']  +  d  <  s,  condition  (a)  implies  that  {/?'’'(!  — e*)}Q  ((7I1  <  m-[A;']  +  d) 
are  bounded  in  £00 (HI*^)-  Thus,  we  obtain  that  the  functions  in  (2.7.2)  are  bounded  in  Li{fl),  for 
some  neighborhood  fl  of  00,  uniformly  in  a  €  Z*^. 

Next,  we  show  that  the  functions  in  (2.7.2)  are  bounded  in  £i(B),  for  some  origin-neighborhood 
B.  Here,  recall  first  that,  on  IR‘*\0,  the  Fourier  transform  of  ip  -  <pa(- +  a)  is 

(2.7.3)  (1-®;)?. 
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Since,  by  (b),  1 — e*  has  a  zero  of  order  m  +  1  at  the  origin,  while,  by  (e),  the  order  of  the  distribution 
0  at  the  origin  is  <  m,  it  follows  that  the  representation  (2.7.3)  for  the  Fourier  transform  of 
0  ~  ^a(‘  +  «)  extends  to  the  entire  IR**. 

Assumptions  (d)  and  (b),  when  combined,  prove  that  is  bounded  around 

the  origin  as  long  as  l/?|i  <  m  -  + 1,  and  has  a  singularity  of  order  ib'  +  |/?|i  -  m  -  1  for  larger  /?. 

This  implies  that,  as  long  as  |/?|i  <  m  -  (*:']  +  <i,  D^(l  -  has  a  singularity  at  the  origin  of 

order  <  d,  and  hence  is  integrable  there.  The  required  uniformity  in  the  Li(B)-norm  follows  from 
the  uniformity  in  assumption  (a).  4 

Under  a  slightly  stronger  assumption  on  the  behaviour  of  4>  around  the  origin,  the  assertions 
of  the  last  theorem  can  be  improved  as  follows. 

Corollary  2.7,4.  Let  A  =  (A(a,^))2z<<xE  coefEcieat  matrix  of  (2.2.1),  and  assume  that 

each  of  the  rows  of  .A  is  in  ^i(H).  Let  e*  be  defined  as  before.  Suppose  the  following  conditions 
bold: 

(a,b,c,d)  Same  as  in  Theorem  2.7.1,  only  that  we  require  m  <  s  -  1. 

Around  the  origin,  <f>  can  be  represented  as  a  sum  po  +  Pi,  such  that 
(do)  On  1R‘^\0,  Po  is  continuous  and  is  a  homogeneous  function  of  order  —k',  i.e.,  for  every  t  >  0, 
the  support  of  the  distribution  Otpa  —  Pa  is  the  origin. 

(di)  Each  D^'Pi,  |q|i  <  r,  (calculated  on  IR‘^\0)  has  a  singular!  ty  of  order  A:' (/Ji)  +  la|i  at  the  origin, 
for  some  positive  k'(pi)  <  [h']. 

(e)  Each  pi  is  a  tempered  distribution  of  order  <  m. 

Then  Condition  2.2.5  is  satisfied  with  mA  ~  m  -  k'  +  d  +  1,  provided  that  mA  <  r,  s. 

The  improvement  here  can  be  observed  as  follows:  assumptions  (do),  (di)  imply  that  our 
function  $  has  a  singularity  of  order  k'  at  the  origin.  Therefore,  a  direct  application  of  Theorem 
2.7.1  can  establish  Condition  2.2.5  only  for  ;=  m  -  [A:']  +  d.  However,  the  present  corollary 
asserts  that  the  Central  Condition  is  valid  also  for  the  larger  tua  m  -  k'  +  d  +  1. 

Proof.  Let  q  be  a  compactly  supported  C°°(]R‘*)-function,  which  is  1  around  the  origin, 
and  such  that  the  decomposition  <^  =  po  +  Pi  is  valid  on  suppij.  We  write 

4>  =  vpq  +  vpi  +  (1  - 

One  easily  checks  that  the  inverse  transform  t  of  qpi  +  (1  —  q)^  satisfies  conditions  (c-e)  of  Theorem 
2.7.1,  with  A:'(r)  =  A:'(pi),  and  hence  with  (A:'(r)J  <  (A:'].  Defining  Tq,  a  G  ZZ*^  by 

To  := 

C€3 

we  may  appeal  to  Theorem  2.7,1  to  conclude  that  Condition  2.2.5  is  satisfied  by  t,  and  with 
i71>i(t)  —  m-  [A:'(pi)]  +d>m  -  [A:']  +  d  +  l>m-fc'  +  d  +  l. 

Consequently,  in  order  to  prove  the  corollary,  we  may  assume  without  loss  that  <l>  =  rjpo,  as 
we  do  in  the  rest  of  the  proof.  This  trivially  implies  that  (f>  satisfies  condition  (c)  of  Theorem  2.7.1 
for  all  r,  (since  is  now  compactly  supported)  and  also  satisfies  (e)  and  (d)  there,  with  k'{<f>)  =  k' 
by  the  assumptions  on  po.  This,  however,  would  lead  us  to  an  unsatisfactory  estimate  for  tua,  and 
thus,  a  finer  analysis  is  required. 

We  let  Tq  be  the  homogeneous  polynomial  of  degree  m  +  1  in  the  Taylor  expansion  of  1  —  e* 
around  the  origin,  and  write 

(1  -  e*  )0  =  (1  -  eDqpo  =  TaVPo  +  (1  -  e*  -  Ta)vpo. 
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Since  po  is  a  distribution  of  order  <  m,  and  is  homogeneous  of  order  -k',  and  since  Tq  is  homo¬ 
geneous  of  order  m  +  1,  the  product  TaPo  is  a  homogeneous  function  of  order  m  +  1  -  k'.  Thus, 
its  inverse  transform  g  decays  at  oo  like  OQ  .  {TaVPaY  *9,  and  decays 

rapidly,  we  conclude  that  {TarjpoY  =  0(1  ■  +<*)). 

The  other  term  in  the  representation  of  (1  -  e*  )(^  is  (1  -  e*  -  Ta)4>.  Here,  all  derivatives  of 
1  -  e*  -  Tq  of  orders  <  m  -f-  1  <  s  vanish  at  the  origin.  Therefore,  the  argument  used  in  the 
proof  of  Theorem  2.7.1  can  be  followed  (with  (1  —  e*)<?i  replaced  by  (1  -  e*  -  Ta)(p)  to  yield  that 
((1  —  e*  —  Ta)4>y  decays  at  00  like  0(|  •  |-("‘+i-t*'l+‘i))^  which  is  (at  worst)  the  required  rate. 

The  uniformity  required  in  Condition  2.2.5  follows,  once  again,  from  the  uniformity  in  condition 

(a)  here.  4^ 

Theorem  2.7.1  as  well  as  Corollary  2.7.4  require  the  matrix  A  to  satisfy  conditions  that  may 
impose  constraints  on  the  distribution  of  E,  and  may  be  hard  to  obtain.  However,  it  is  shown  in 
[BuDL]  (see  also  §3.4)  that,  for  every  fi.xed  m,  and  for  all  sufficiently  small  c  <  co(m)  if  each  ball  of 
radius  c  in  intersects  with  E,  there  exists  a  matrix  A  that  satisfies  conditions  (a-b)  of  Theorem 
2.7.1. 


2.8.  On  the  decay  rates  m^, 

The  derivation  of  the  main  results  of  this  section  required  us  to  assume  decay  rates  on  the 
locali2ation  sequence  p,  the  localization  function  ip,  and  the  difference  4>{-  —  a)  —  <pa  (cf.  Conditions 
2.1.2  and  2.2.5).  We  preferred  to  write  these  assumptions  in  terms  of  the  0(|  |~*^)  (for  an  appropriate 
k)  since  that  allows  us  to  derive  quantitative  results  concerning  the  behaviour  of  {ipci)a&7i*  at  00  (cf. 
Lemma  2.2.6),  which  is  important  in  several  places,  especially  in  §2.4.  However,  as  Theorem  1.2.5 
actually  asserts.  Theorems  2.2.9,  2.2.16,  and  2.3.1  remain  valid  under  slower  decay  rates.  These 
improvements  are  not  very  important  with  respect  to  the  decay  rates  of  p  (for  example,  the  decay 
rates  on  p  assumed  in  Theorem  2.2.9  are  more  restrictive  than  those  assumed  in  the  Localization 
Conditions  2.1.2).  Also,  such  improvements  do  not  seem  to  be  very  significant  with  respect  to  the 
decay  rates  of  ip  (in  2.1.2),  since,  in  the  known  examples,  ip  decays  at  00  at  faster  rates.  However, 
as  far  eis  the  behaviour  of  <p{-  -a)-<pc  is  concerned,  every  relaxation  of  the  Central  Condition  2.2.5 
is  important,  since  it  might  simplify  the  construction  of  the  matrix  A,  hence  eventually  simplify  the 
actual  approximation  scheme  La-  Therefore,  we  devote  the  present  subsection  to  a  brief  discussion 
of  these  relaxations. 

Our  modified  conditions  are  as  follows: 

(a)  In  lieu  of  (a)  in  Conditions  2.1.2,  we  assume  only  that  p  € 

(b)  In  lieu  of  (b)  in  Conditions  2.1.2,  we  assume  only  that  “  “)l  ^  ^oo(R‘^)- 

(c)  In  lieu  of  Condition  2.2.5,  we  assume  that  ~  ~  M  ^  ■too(R‘*)- 

The  new  variant  of  Lemma  2.2.6  then  reads  as  follows: 

Lemma  2.2.6*.  Assume  that  conditions  (a-c)  above  bold.  Then 

X;  IV'al  €  Ioo(R‘'), 
a€2Z'' 


and,  hence,  the  operator 

c-*  i’acia) 

a€Z« 
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is  bounded  from  tooi'S,'^)  to  LoollR-*^)- 

Proof.  By  the  definition  of  xj;  and  xpa, 

iV'aW  -  rlj{x  -  a)l  <  |/t()9  -  a)\\<f>^{x)  -  ^(i  -  0)\. 

^€22-* 


Therefore, 

52  Iti'„(x)  -  Xp{x  ~a)\<  52  -  4>{x  -  0)\  <  llMllr,(2-«)  5Z  ~ 

a€Z2''  a,0€Z'‘  SeZ"* 

Consequently,  by  (a)  and  (c)  above,  ~  ^)1  ^  and  combining  that  with 

(b)  above  we  obtain  the  lemma’s  claim.  4 

The  results  of  Theorems  2.2.9,  2.2.16,  and  2.3.1  remain  unchanged,  and  with  the  same  proofs. 
In  fact,  one  observes  that  Condition  2.2.5  is  not  invoked  in  the  proof  of  Theorem  2.2.9  in  its  full 
power,  and  only  the  Loo-boundedness  of  Yla^Ti*  ~  ^  there. 


3.  Examples 

3.1.  The  approximation  schemes  of  [BR] 

The  reference  [BR]  contains  a  general  study  of  Loo(ni‘^)-approximation  orders  from  spaces 
generated  by  the  shifts  of  one  basis  functions  It  assumes  xp  to  decay  fast  enough  at  oo  to  make 
the  map 

xp*'  :c^xp**  c:=  52  “  0()c(a) 

aez" 

well-defined  and  continuous  from  to  ZrooClR*^)  (see  (1.2.2)).  Note  that  such  a  condition 

implies  that  xp  G  I’i(IR‘^),  and  is  implied  by  (b)  of  Localization  Conditions  2.1.2.  After  establishing 
upper  bounds  on  approximation  orders  from  the  “span”  of  the  shifts  of  xp  (cf.  Theorem  2.8  there,  and 
also  Lemma  2.2.15  above),  [BR]  attempts  to  realize  these  bounds  by  employing  the  appro.ximation 
map 

(3.1.1)  r.:f^xP*'A(f)\^^, 
where  A  is  the  convolution  map  A  ♦  /,  with 

(3.1.2)  A  =  ri/xp. 

Here,  xj  is  any  smooth  compactly  supported  function  which  is  1  around  the  origin.  Results  which 
are  relevant  to  the  basis  functions  of  this  paper  are  obtained  in  Theorems  3.6  and  4.2  of  [BR].  In 
particular,  the  following  result  follows  from  Theorem  4.2  there: 

Result  3.1.3.  Assume  that  xp*'  is  bounded,  xp{0)  ^  0,  and  that,  on  IR‘^\0,  xp  =  pixp,  with  ft 
continuous  and  2ir-periodic  and  (p  continuous  on  IR‘^\0.  Assume  further  that  for  some  neighborhood 
Q  of  the  origin  and  some  positive  k, 

52  ll^('  +  °f)lli.oo(n)  < 

a62ir2Z''\0 
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8Lnd 

III  •  r*/^lkoo(n)  <  <x). 

Then  L  provides  approximation  order  (at  least)  k  for  every  /  e  U^^(1R‘*). 

In  order  to  extend  Result  3.1.3  to  non-uniform  grids,  we  want  to  apply  Theorem  2.3.1,  hence 
need  to  verify  its  various  ronditions.  First,  A  of  the  present  L  is  indeed  a  convolution  operator  A*. 
Since  the  result  assumes  V'  to  be  non-zero  at  0,  and  -0  is  continuous  (since  rp  G  Li(lR‘*)  by  virtue 
of  the  boundedness  of  ^*'),  we  conclude  that  A  is  bounded.  Furthermore,  around  the  origin 

1  |x)"*';2(x))  =  ){|ir'=/^(x))'0(x)l  <  const, 

and  hence  |  ■  is  bounded  (on  IR*^).  Therefore,  among  the  various  conditions  assumed  for  (a)  of 
Theorem  2.3.1,  we  need  only  assume  the  Localization  Conditions  2.1.2  and  the  Centred  Condition 
2.2.5; 

Corollary  3.1.4.  Let  L  be  defined  as  in  (3.1.1)  and  (3.1.2),  and  assume  that  ^  satisfies  the 
assumptions  of  Result  3.1.3.  If,  in  addition,  the  Conditions  2.1.2  an^2.2.5  are  fulfilled,  the  corre¬ 
sponding  La  provides  approximation  order  at  least  k  to  every  f  e  1'F^(1R‘^),  with  k  as  in  Result 
3.1.3. 

Under  the  further  assumption  that  5^  0,  [BR]  proves  that  (3.1.1)  attains  the  highest 

possible  approximation  order.  Thus,  the  above  corollary  shows  that,  under  all  the  assumptions 
made  there,  we  are  able  to  realize  in  the  scattered  case  the  best  possible  approximation  orders 
available  in  the  uniform  case. 


3.2.  Interpolation  Schemes 

Approximation  via  cardinal  interpolation  that  employs  radial  basis  functions  is  studied  in 
Buhmann’s  thesis  (cf.  [Bul,2|).  The  interpolation  operator  is  of  the  form 

Lf=^  rp(-~a)f(a), 
a€2Z<‘ 

with  Ip  the  fundamental  solution  of  the  cardinal  interpolation  problem,  i.e.,  a  function  that  van¬ 
ishes  on  7L\Q,  and  assumes  the  value  1  at  the  origin.  The  function  tp  is  obtained  as  an  infinite 
combination  of  shifts  of  the  original  basis  function  (p: 

The  sequence  p,  though  not  constructed  for  the  mere  sake  of  localization,  is  indeed  a  localization 
sequence  since  tp  is  shown  to  decay  at  oo,  in  contrast  with  the  original  <p. 

The  existence  of  the  above  fundamental  solution  is  proved  in  (Bul,l^  under  certain  condi¬ 
tions  on  the  Fourier  transform  <p  of  the  basis  function  <p.  In  particular,  <p  is  assumed  to  be  ra¬ 
dially  symmetric,  to  be  smooth  and  positive  on  ]R**\0,  to  have  derivatives  that  decay  at  oo  at  a 
0(1  •|~f'''''^l)-rate,  and,  most  importantly,  to  have  a  singularity  of  a  certain  type  and  of  some  positive 
order  k'  at  the  origin.  The  sequence  p  is  then  proved  there  to  decay  at  oo  at  a  rate  0(1  •  )), 

and  the  same  decay  rate  is  then  established  with  respect  to  the  fundamental  interpolant  ip.  The 
sum  that  defines  ip  is  proved  to  converge  uniformly  on  compact  sets.  Thus,  we  see  that  the  Local¬ 
ization  Conditions  2.1.2  as  well  as  the  Uniform  Scheme  Condition  2.1.6  which  are  required  in  our 
conversion  theorems,  are  valid  in  Buhmann’s  interpolation  schemes. 
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The  approximation  orders  associated  with  cardinal  interpolation  that  are  derived  in  [Bui, 2] 
are  of  the  following  type.  In  case  k'  =  k  +  e,  with  0  <  e  <  \,  k  integer,  L  is  proved  to  provide 
approximation  order  k'  to  functions  in  C*'''^(IR‘^)  whose  derivatives  of  order  A:,  it  +  1  are  bounded. 
For  an  integer  k',  L  is  proved  to  satisfy  IWh/ -L{ahf)\\ao  =  log  /i|),  for  functions  in 

with  bounded  derivatives  of  order  k'  -  l,Jt'.  Furthermore,  for  an  even  integer  k',  and  under  a 
stronger  assumption  on  the  behaviour  of  <f)  at  the  origin,  the  0(/i*  |  log /i|)  order  is  improved  to 
approximation  order  k'.  Though  we  do  not  provide  here  the  details  of  those  extra  conditions  (nor 
we  provide  the  full  details  of  Buhmann’s  original  conditions),  we  mention  that,  with  the  aid  of  the 
extra  conditions,  is  proved  to  decay  like  0(1  •  for  some  positive  integer  u. 

In  the  conversion  of  the  above  results  to  the  scattered  case,  we  need  not  appeal  to  the  results  of 
section  2.3,  simply  because  /i  decays  here  at  a  rate  which  allows  the  application  of  Theorem  2.2.16 
in  its  full  power.  Inched,  we  already  realized  that  Conditions  2.1.2  and  2.1.6  are  valid  here,  and 
the  requirements  on  <f)  needed  in  [Bui, 2]  are  by  far  more  restrictive  than  our  mere  assumption  (a) 
of  Theorem  2.2.16.  The  decay  conditions  on  fi  required  in  (c)  of  Theorem  2.2.16  are  satisfied  only 
in  the  case  of  a  non-integer  k'.  Still,  this  is  in  agreement  with  Buhmann’s  results,  since  he  obtains 
only  the  order  0(/i*''  "^^1  log/i|)  in  the  integer  case,  exactly  as  can  be  derived  for  (cf.  the  remark 
after  the  proof  of  Theorem  2.2.16).  Finally,  in  the  case  when  k'  is  an  integer  and,  nonetheless,  a 
full  approximation  order  is  obtained  in  [Bul,2],  ft  decays  fast  enough  to  be  well-defined  on  11*;', 
and  Theorem  2.2.16  applies  again  here.  In  summary,  assuming  that  the  Central  Condition  2.2.5  is 
satisfied,  all  the  approximation  orders  that  arc  established  in  [Bui, 2]  for  cardinal  interpolation  L 
are  valid  also  for  its  “scattered  variant”  (though,  L^,  of  course,  is  not  a  cardinal  interpolant, 
nor  is  guaranteed  to  be  any  kind  of  interpolant),  when  the  approximand  satisfies  the  conditions 
required  in  [Bui, 2],  and,  in  addition,  is  bounded. 

An  extension  of  the  £/i-approximation  orders  to  unbounded  functions,  follows  from  the  analysis 
of  §2.4.  For  example,  if  k'  =  k  +  e,  Q  <  £  <  1,  then,  the  decay  rate  of  and  of  /r  are  = 

d  +  k',  which  allows  us  to  invoke  Theorem  2.4.4  for  the  value  r  :=  k  there,  provided  that  of 
Condition  2.2.5  is  known  to  be  larger  than  d  +  k.  Consequently,  the  approximation  order  k'  of  La 
holds,  by  Theorem  2.4.4,  not  only  for  bounded  functions,  but  for  functions  in  C*'''^(]R‘^)  whose 
derivatives  of  orders  k  and  fc  -I-  1  are  bounded.  As  mentioned  before,  this  is  exactly  the  class  of 
functions  that  arc  approximated  in  the  uniform  grid  results  of  [Bui, 2].  Thus,  the  conversion  here 
fully  preserves  both  the  approximation  order  and  the  space  of  approximands.  Analogous  results 
are  valid  in  the  case  of  an  integer  k'. 

3.3.  Quasi-interpolation  schemes 

Quasi-interpolation  schemes  which  are  based  on  uniform  shifts  of  localized  radial  basis  functions 
are  discussed  in  [J],  [DJLR],  [P],  [Ra],  [Bui, 3].  In  this  section,  we  show  that  the  uniform  grid  results 
of  [DJLR]  can  be  extended  to  non-uniform  grids,  with  the  aid  of  the  machinery  of  §2  here. 

In  [DJLR],  the  Fourier  transform  of  the  basis  function  <p  is  assumed  to  be  represented,  on 
1R'^\0,  as  a  quotient 

(3.3.1)  ^  =  F/G, 

with  F,  G  satisfying  several  conditions,  among  which  we  mention  here  the  following  two: 

(a)  G  is  a  homogeneous  polynomial  of  degree  A:'  e  2  IN,  with  no  zeros  in  1R‘*\0. 

(b)  F  €  G°°(1R‘^\0),  is  smooth  (but  not  necessarily  infinitely  smooth)  at  the  origin  and  does  not 
vanish  there. 
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We  forgo  mentioning  the  detailed  conditions  assumed  in  [DJLR]  (and  refer  the  reader  to  that 
reference),  since  most  of  those  details  will  not  be  needed  here.  We  do  mention  that  the  basic 
examples  of  <f)  that  satisfy  all  the  conditions  assumed  there  are  the  fundamental  solutions  of  the 
iterated  Laplacian  (in  particular,  the  thin-plate  splines  in  even  dimensions)  and  their  “shifted” 
version,  where  |a:|  is  replaced  by  (jip  +  c^)^/2  (in  particular,  the  multiquadrics  in  odd  dimensions). 

The  localization  process  of  [DJLR]  is  done  with  the  aid  of  a  finitely  supported  localization 
sequence  fi.  This  means  that,  for  the  sake  of  extending  the  [DJLRj-schemes  to  the  scattered  case, 
the  required  conditions  (a)  and  (c)  of  2.1.2  are  automatically  satisfied.  Since  the  error  analysis  of 
[DJLR]  is  based  on  polynomial  reproduction,  a  careful  attention  is  given  there  to  the  decay  rates 
of  Ip  at  oo.  In  general,  these  decay  rates  are  =  d  +  f  +  1  for  some  non-negative  integer  i, 
which  means  that  condition  (b)  of  2.1.2  is  also  satisfied  in  the  present  case. 

The  approximation  scheme  used  in  [DJLR]  is  the  simple  one 

(3.3.2)  L/;=  Y. 

and  the  quasi-interpolation  argument  is  employed  there  to  yield  the  error  estimate 

(3.3.3)  \M  -  Lcr^fW  =  0(/i^'+*l  log  h\),  (!  :=  min{£,  *'-!}, 

provided  that  /  is  ^'-M-times  differentiable,  and  that  |/|f',oo  and  |/|f'+i,oo  are  finite.  For  functions 
/  whose  derivatives  of  orders  £'4-2  are  also  continuous  and  bounded,  and  under  more  subtle 
information  on  the  decay  of  tp,  the  above  0{h^  ’♦■^llog/il)  is  improved  to  0{h^''^^).  We  refer  to 
section  4  of  [DJLR]  for  the  full  details. 

It  is  now  easy  to  verify  that  all  the  conditions  on  (p,  rp,  and  p  imposed  in  Theorem  2.2.9  and 
Theorem  2.2.16  are  satisfied  by  the  schemes  considered  in  [DJLR].  Therefore,  the  only  condition 
that  we  really  need  to  assume  is  the  Central  Condition  2.2.5,  in  order  to  obtain 

Theorem  3.3.4.  Let  <p  be  any  of  the  basis  functions  considered  in  [DJLR]  (whose  Fourier  trans¬ 
form  has  a  singularity  of  even  order  k'  at  the  origin).  Let  ip  be  its  localization  (localized  with  the 
aid  of  a  finite  p,),  which  decays  at  oo  like  C?(|  •  £,et  £,  be  the  approximation  scheme  of 

(3.3.2),  and  let  La  be  its  scattered  version.  If  Condition  2.2.5  is  satisfied,  then 

(((L  -  =  0(A"), 


for  every  f  €  W,^(IR'^)  D  C*'(]R‘^),  provided  that  k  <  k'. 

Note  that  the  result  is  valid  (or  k  :=  k',  even  when  t  <k'. 

Thus,  while  the  approximation  orders  of  (3.3.3)  and  their  variants  are  restricted  by  the  rate 
of  decay  of  tp,  i.e.,  the  parameter  I,  the  estimate  for  (L  —  LA){(rhf)  is  optimal  already  for  £  =  0. 
The  extension  in  the  above  theorem  from  uniform  grids  to  scattered  grids  required  us  to  pay  some 
price:  while  the  uniform  results  are  obtained  with  respect  to  functions  whose  derivatives  of  order  £ 
and  £4-1  are  continuous  and  bounded,  the  non-uniform  extensions  require  /  to  be  bounded,  too. 
However,  this  further  restriction  can  be  relaxed  with  the  aid  of  Theorem  2.4.4:  choosing  k  and  r  in 
that  theorem  to  be  our  £  here,  we  see  that  since  in  our  present  case  =  oo,  and  =  d  4-  £  4- 1, 
we  need  only  know  that  in  Central  Condition  2.2.5  ruA  >  d  4-  £,  in  order  to  obtain  that 

for  functions  /  whose  derivatives  of  order  £  and  £  -H  are  continuous  and  bounded. 


25 


Our  last  remark  here  deals  with  the  polynomials  in  Ss(<p)-  It  is  shown  in  [DJLR]  (cf.  §3  there) 
that  if  F  of  (3.3.1)  is  the  constant  function  (as  is  the  case  for  the  thin-plate  splines),  then  we  can 
construct  the  function  V’  to  have  an  arbitrary  high  decay  rate  and  to  reproduce  n^fTP,  with  V 
the  kernel  of  G{D),  and  k  =  -  d  -  1.  It  is  further  shown  there  that  the  corresponding  sequence 

fi  would  then  annihilate  11*  n  P.  In  such  a  case,  Corollary  2.5.3  can  be  invoked  to  show  that  La 
reproduces  11*  flP,  and,  consequently,  that  contains  11  DP. 

3.4.  Review  and  extensions  of  the  scheme  used  in  [BuDL] 

As  we  already  mentioned  in  the  introduction,  the  reference  [BuDL]  discusses  approximation 
schemes  and  approximation  orders  that  are  based  on  the  scattered  shifts  of  one  beisic  function  <f>. 
In  this  subsection  we  review  the  results  of  [BuDL]  with  the  aid  of  the  techniques  and  observations 
of  the  present  paper  and,  in  this  course,  extend  those  results  in  several  directions. 

The  aim  of  [BuDL]  is  to  extend  the  approximation  orders  on  uniform  grids  of  [DJLR]  to  non- 
uniform  grids.  While  we  did  the  same  in  the  previous  subsection,  the  error  analysis  we  have  used 
here  and  the  one  used  in  [BuDL]  are  entirely  different.  Furthermore,  apparently,  the  approximation 
schemes  used  here  and  there  seem  to  be  very  different,  too.  We  will  show,  however,  that  the  [BuDL] 
scheme  can  be  viewed  as  a  variation  of  a  special  case  of  the  method  we  suggest  here.  In  particular, 
implicitly,  the  approximation  scheme  of  [BuDL]  is  engaged,  too,  with  the  approximation  of  uniform 
shifts  by  scattered  shifts. 

The  setup  in  [BuDL]  starts  with  a  basis  function  (f)  that  obeys  the  conditions  of  [DJLR]  (some 
of  these  conditions  are  quoted  in  the  previous  subsection).  Then,  for  each  center  ^  6  H,  a  function 
of  the  form 

»?e2 

with  (/i^, ,,),,€=  a  sequence  of  finite  support  (supported  on  centers  which  are  close  to  ^),  is  con¬ 
structed.  The  functions  (V’$)«e3  are  constructed  such  that  they  decay  (uniformly)  at  oo  at  rates 
similar  to  those  of  V'  of  [DJLR],  i.e.,  0(1  •  j-‘'-*“i)  for  some  non-negative  integer  t.  The  approxi¬ 
mation  map  La  is  then  defined  as 


A  very  careful  choice  of  the  coefficients  (m{,ij)  enables  [BuDL]  to  extend  the  results  of  [DJLR]  to 
the  scattered  case. 

We  make  two  preliminary  observations  coiwerning  the  difference  between  the  map  La  of 
[BuDL]  and  the  map  La  here.  First,  the  map  La  uses  the  values  of  /  at  E  only,  while  La  in¬ 
volves  evaluation  at  71^  of  A  *  /.  Second,  the  localization  sequence  which  is  used  in  [DJLR]  in 
order  to  construct  the  localization  V’  from  the  original  <!>  (and  which  is  used  in  the  approach  of  the 
present  paper  in  an  essential  way)  is  not  invoked  directly  in  [BuDL].  Instead,  it  is  imitated  in  the 
sequences  (/i{,„)„€E,  ^  €  E. 

The  problem  of  selecting  appropriate  coefficients  (/ij.jj)  is  approached  in  [BuDL]  as  follows. 
First,  it  assumes  the  existence  of  a  matrix 

~  (-^(®>C))a€K^,€e3 

that  satisfies  the  following  three  conditions: 
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Assumptions  3.4.1. 

(a)  There  exists  a  positive  constant  ci  such  that  A(a,  0  =  0,  whenever  |a  -  ^1  >  cj . 

(b)  The  sequence  set  :  a  €  E*'},  lies  in  /i(E)  and  is  bounded  there. 

(c)  For  some  positive  integer  m,  and  for  every  a  G  2Z‘',  h  ~  n^,  i  e., 

E«e3  =  P{o),  for  every  pe  !!„. 

Note  that  the  above  assumptions  imply  conditions  (a-b)  of  Theorem  2.7.1.  Indeed,  most  of  the 
subsequent  results  here  that  assume  Assumptions  3.4.1,  are  valid  under  weaker  assumptions  (such 
as  (a-b)  of  Theorem  2.7.1). 

The  value  of  m  required  in  [BuDL]  is  m  =  A:'  4-^,  with  i  as  above,  and  k'  the  degree  of  G  in  (3.3.1). 
The  coefficients  M  :=  (m^.ij)  are  then  defined  by 

Aia,aA(p,v){N{--a).N{--0)), 

a,^€7L* 


where  (•,  •)  is  a  certain  bilinear-form  that  depends  on  4>,  and  where  iV  is  a  sufficiently  smooth 
compactly  supported  multivariate  spline  which  reproduces  H^'+r  (i.e.,  Yae7i<^  ~  ^)p(“)  ~  Pi 
Vp  €  11^'+/),  (cf.  [BuDL]  for  the  full  details).  With  M  being  the  Gramian  matrix 

(3.4.2)  M  :=  {uia  -  /3)),,p  :=  {{N{-  -  a),  N{-  -  /?))), 


we  see  that 

M  =  A'^AfA. 


Furthermore,  it  is  proved  in  [BuDLj  that  the  sequence  i/  is  a  suitable  localization  sequence  for  the 
[DJLR]  construction. 

The  crucial  properties  of  M  required  in  [BuDL]  are: 

(i)  Ynes  Pc.v  P(^)  ~  every  p  G  fl  V,  and  every  ^  G  S,  where  V  is  the  kernel  of  the 
differential  operator  G{D),  with  G  the  denominator  of  (f>  (see  (3.3.1)). 

(“)  E4€hP«.»/P(0  =  0i  for  every  p  G  and  every  p  G  E. 

The  [BuDL]  analysis  insists  on  a  symmetric  p,  since,  in  that  way,  condition  (ii)  above  is  implied 
by  (i),  eis  G  is  a  homogeneous  polynomial  of  degree  k'. 

From  our  factorization  A^AfA  of  M,  one  can  observe  that  (i)  and  (ii)  hold  with  A/"  =  ip{a  — P))a,0i 
for  any  localization  sequence  p  that  satisfies  the  [DJLR]  requirements,  and  not  just  for  the  specific 
}/,  (3.4.2),  that  is  employed  by  [BuDL]. 

In  order  to  see  how  the  [BuDL]-scheme  is  related  to  the  approach  of  the  present  paper,  we 
denote 

^  :=  V’E  =  ^  :=  h(<^)  =  {<!>{•  -0-  ^  e  E}. 


Then,  we  may  write 


’J'  = 


Also,  since  v  can  serve  as  a  localization  sequence  of  [DJLR],  it  follows  that 


^:=  ^  t/(Q)<^(--Q) 

aeZ'' 


decays  at  oo  like  0(|  •  Therefore,  we  see  that  the  Localization  Conditions  2.1.2  are 

satisfied  for  the  choice  p  :=  i/. 
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Given,  now,  an  approximand  /,  we  may  write  the  approximant  La/  that  is  provided  by  [BuDL] 
in  matrix  form  as 

with  /|.  the  restriction  of  /  to  E,  treated  as  a  column  vector.  On  the  other  hand,  the  approximation 
map  La  of  the  present  paper,  when  written  in  matrix  form,  too,  reads,  in  case  the  convolution  A* 
is  the  identity,  as 

LAf  =  (f\^yA/’A^, 

since  are  exactly  our  pseudo-shifts  {<t>a)ai  and  AfA^  are  our  pseudo-shifts  (^a)a- 
In  summary,  while  our  map  La  is  of  the  form 

La!  = 

ae2Z‘' 

the  [BuDL]  map,  written  with  the  aid  of  our  {‘>pa),  is  of  the  form 

LaI  =  (Af^^){a)'ipa- 

Consequently,  as  soon  as  we  can  show  that  provides  a  “reasonable”  approximation  to  , 
we  will  be  able  to  derive  the  [BuDL]  approximation  results  directly  from  the  results  of  the  present 
paper.  Since  the  results  of  the  present  paper  are  derived  under  assumptions  which  are  substantially 
more  general  than  the  ones  of  [BuDL],  we  will  obtain  in  such  a  way  an  extension  of  the  ideas  and 
the  results  of  that  paper  beyond  their  stated  limits. 

The  fact  that  (/l/|,)(a)  approximates  /(a)  and  in  the  “right”  way,  can  be  established  in  a  far 
more  general  setup  than  the  one  discussed  above,  as  we  now  describe.  First,  we  adopt  Assumptions 
3.4.1,  with  a  general  integer  m  whose  desired  values  will  be  clarified  in  the  sequel.  Next,  let 
/  €  and  a  €  2'^  Let  Taf  be  the  Taylor  polynomial  of  degree  m  of  /  about  a.  By 

Assumption  3.4.1  (c),  A(p|2)  =  for  every  p  €  Dm,  and  therefore 

A((T,/),3)(a)  =  (To/)(a)  =  /(a). 

This  implies  that,  in  view  of  Assumptions  3.4.1  (a,b)  on  A, 

iWi.Ko)-/(a)i=i  Y, 

\a-il<ci 

^COnSt„(^ci  1|A(Q!,  (2j|/(oo,m+l  —  const  |/|oo,m+l 

Consequently,  we  obtain 

Lemma  3.4.3.  For  /  e  under  Assumptions  3.4.1,  Localization  Conditions  2.1.2,  and 

Central  Condition  2.2.5, 

HLa  ~  L/l)((Tft/)|  <  const /l”*'^^l/|oo,m+l • 

Proof.  From  the  arguments  that  precede  the  lemma,  we  conclude  that 

1(£<4  —  L/i)(<T/i/)I  <  const /l”*'^^|/|oo,m+l  ^  ]  |'0a(l)|  <  const  A  I/loo, m+lj 

oez'* 

where,  in  the  second  inequality,  we  used  the  fact  that  IV'al  €  Zroo(IR-‘*)i  which,  thanks  to 

our  assumptions  here,  follows  from  Lemma  2.2.6.  ♦ 

Therefore,  in  order  to  extend  the  approximation  order  results  concerning  the  map  La  to  the 
map  La  we  need  the  number  m+1  of  Lemma  3.4.3  to  match  or  exceed  the  orders  for  ||(Xyi— Zr)(«7/i/)|| 
derived  in  the  theorems  of  section  2.  This^can  be  guaranteed  if  we  take  m  as  the  largest  integer 
smaller  than  the  order  of  the  singularity^  ^  has  at  the  origin.  For  example,  the  following  is  the 
extension  of  Theorem  2.2.16  to  the  map  La'- 
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Corollary  3.4.4.  Assume  that  conditions  (a~d)  of  Theorem  2.2.16  hold.  Assume  further  that  the 
matrix  A  satisfies  Assumptions  3.4.1,  with  m  taken  to  be  k  of  Theorem  2.2.16.  Then,  the  assertions 
of  that  theorem  remain  valid  with  respect  to  the  map  defined  above. 

Proof.  Since  all  the  approximation  orders  asserted  in  Theorem  2.2.16  with  respect  to  La 
are  <  A;  +  1,  the  present  assertion  follows  from  Lemma  3.4.3. 

Remark.  With  k,k'  as  in  Theorem  2.2.16,  if  k'  ^  ZZ+,  then  Condition  2.2.5  (that  is,  condition 
(d)  in  Theorem  2.2.16),  which  is  the  crucial  condition  exploited  in  §2,  becomes  redundant  in  the 
above  corollary,  under  mild  conditions  on  4>.  This  follows  by  an  application  of  Corollary  2.7.4: 
conditions  (a)  and  (b)  of  that  corollary  are  a  consequence  of  Assumptions  3.4.1,  and  therefore,  if  <j) 
satisfies  assumptions  (c-e)  of  that  corollary,  we  obtain  that  of  Condition  2.2.5  can  be  taken  as 
m  —  k'  +  d  +  l  =  k  —  k'  +  d  +  l>d,  since  k  —  k'  >  —1,  as  fc'  ^  7Z+.  In  the  other  case,  when  k'  is 
an  integer,  we  obtain  a  similar  result  if  m  of  Assumption  (c)  of  3.4.1  is  taken  as  k'. 

We  show  now,  that  the  approximation  order  results  of  [BuDL]  follow  from  our  results  in  §3.3 
and  §3.4.  Indeed,  for  bounded  approximands  this  is  obtained  directly  from  Theorem  3.3.4  and 
Corollary  3.4.4.  As  for  unbounded  approximands,  the  reasoning  goes  as  follows.  First,  by  the 
discussion  in  the  last  paragraph  of  §3.3,  our  scattered  analog  La  of  the  [DJLR]  L  provides  to 
functions  with  bounded  f-  and  I  +  1-order  derivatives  the  desired  approximation  orders,  provided 
that  niA  of  Condition  2.2.5  is  >  d  -f-  f .  In  particular,  La  is  well-defined  on  Ni,  and,  adopting  (a) 
and  (b)  of  Assumptions  3.4.1,  we  easily  conclude  that  La  is  well-defined  on  Nt,  as  well.  Thus, 
by  Lemma  3.4.3,  La  provides  the  desired  approximation  orders  as  long  as  m  is  at  least  £.  In 
summary,  the  [BuDL]-scheme  approximates  unbounded  functions  to  the  desired  rates,  as  soon  as 
(i):  tua  >  d  +  i,  (ii)  m  of  Assumptions  3.4.1  is  at  least  i. 

We  can  obtain  (i)  above  from  Assumptions  3.4.1  by  invoking  Corollary  2.7.4,  which  applies  to 
ail  <p  considered  in  [DJLR]  for  the  choice  po  :=  F{^)IG  (with  F,G  as  in  (3.3.1)).  Condition  (i) 
follows  from  that  corollary  ifm  —  fc'-hd-l-l  >  d-\-  £,  i.e.,  if  m  >  fc'  -f-  f.  Such  condition  covers  the 
required  (ii),  as  well.  Hence,  the  choice  m  :=  A:' -I- f  in  Assumptions  3.4.1  (which  indeed  is  the  value 
of  m  used  in  [BuDL])  suffices  to  reproduce  the  [BuDL]  results. 

We  have  just  shown  that  Corollary  3.4.4  reproduces  the  results  of  [BuDL].  In  fact,  it  extends 
those  results  quite  substantially.  First,  it  applies  to  a  much  larger  family  of  basis  functions  <f>,  and 
in  particular,  the  order  of  the  singularity  of  cj)  need  not  be  an  even  integer.  Second,  neither  do  we 
need  the  localization  sequence  p  to  be  the  specific  one  employed  in  [BuDL],  nor  we  even  need  p 
to  be  finitely  supported.  Third,  though  we  did  not  take  advantage  of  that  option,  the  matrix  A 
that  is  used  to  define  i<l>a)a  iiccd  not  be  identical  to  the  matrix  A  used  in  the  expression  (A/(-)(a) 
(they  only  need  to  satisfy,  both.  Assumptions  3.4.1  with  m  =  k'  + 1,  and  m  =  I  respectively).  Also, 
for  bounded  approximands,  requirement  (c)  in  Assumptions  3.4.1  is  weaker  in  the  last  corollary,  as 
compared  to  the  analogous  requirement  in  [BuDL]  (we  need  m  —  I  instead  oivti  —  k'  + 1). 
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